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CHAPTER  I 
INTRODUCTION 


An  adaptive  array  is  an  antenna  system  that  can  automatically 
adjust  its  antenna  pattern.  The  antenna  pattern  is  adjusted  to  track 
desired  signals  and  at  the  same  time,  to  steer  nulls  toward  interfer¬ 
ing  signals.  Such  antennas  are  especially  useful  for  protecting 
communication  links  from  interference  or  jamming,  and  have  both  com¬ 
mercial  and  military  applications. 

Several  techniques  for  controlling  the  weights  in  an  adaptive 
array  have  been  studied.  These  include  the  LMS  array[l],  the  Apple- 
baum  array[2],  and  the  Shor  array[3]. 

This  thesis  is  concerned  with  the  use  of  the  LMS  adaptive  array 
in  a  frequency-shift  keyed  (FSK)  communication  system[4].  The  LMS 
array  has  previously  been  used  in  several  types  of  communication 
systems.  For  example,  the  LMS  array  was  used  in  an  early  adaptive 
system  using  amplitude  moduation[5] .  In  another  experimental  system, 
binary  information  was  biphase  modulated  and  received  by  an  LMS 
array[6]. 

The  communication  system  that  we  will  examine  is  a  binary  fre¬ 
quency  shift  keyed  system[4].  We  assume  a  sequency  of  bits,  l's  and 
0's,  is  to  be  transmitted.  As  illustrated  in  Figure  I  - 1 ,  a  0  is 
transmitted  by  sending  a  burst  of  signal  at  frequency  w0,  and  a  1  by 
a  burst  signal  at  frequency  wj .  In  either  case,  the  time  duration  of 
a  bit  is  T.  Thus  the  transmitted  signal  shifts  back  and  forth  between 
the  two  frequencies,  <d0  and  w| .  At  the  receiver,  bit  detection  is 
done  by  determining,  during  each  bit  interval,  whether  the  received 
signal  has  frequency  w0  or  w| . 

The  adaptive  array  that  we  will  consider  is  the  LMS  array.  The 
LMS  array  adjusts  its  antenna  pattern  to  achieve  least  mean-square 
error  between  the  received  signal  and  a  reference  signal.  Ideally, 
the  reference  signal  should  be  an  exact  replica  of  the  desired  signal. 
However,  in  an  actual  communication  system,  the  received  signal  is 
not  known  in  advance.  Therefore,  the  reference  signal  is  obtained 
by  estimating  the  desired  signal  from  the  array  output. 

The  FSK  system  studied  in  this  thesis  uses  bit  predictions  to 
generate  a  reference  signal.  In  Chapter  II,  a  coding  procedure  is 
presented  which  allows  the  desired  signal  bits  to  be  predicted  at 
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Figure  1-1.  FSK  modulation. 


the  receiver.  Each  bit  is  predicted  from  the  previous  received  bit. 

A  processing  loop,  described  in  Chapter  IV,  uses  the  bit  predictions 
to  process  the  array  output  and  to  form  a  reference  signal.  A  com¬ 
puter  simulation  program  has  been  used  to  demonstrate  the  performance 
of  the  entire  system.  Finally,  Chapter  V  considers  some  possible 
problems  associated  with  the  FSK  system. 


CHAPTER  II 

THE  LMS  ADAPTIVE  ARRAY 


A.  Introduction 

In  this  chapter,  we  will  discuss  the  LMS  algorithm  used  in  many 
adaptive  arrays.  This  discussion  will  provide  some  background  infor¬ 
mation  on  adaptive  arrays  and  also  establish  notation  to  be  used 
later.  We  will  begin  with  the  configuration  of  the  adaptive  array. 
Then  we  will  derive  the  LMS  algorithm  used  to  control  the  array 
weights.  In  the  next  part  of  this  chapter,  we  will  introduce  an 
analytic  signal  model  of  the  LMS  array.  The  use  of  analytic  signals 
makes  the  notation  more  compact  and  easier  to  work  with.  We  will  use 
the  analytic  signal  notation  in  the  following  chapters  to  study  an 
adaptive  array  in  an  FSK  communication  system. 


B.  The  LMS  Algorithm 

A  diagram  of  the  LMS  adaptive  array  is  shewn  in  Figure  1 1 - 1 . 
The  signal  received  by  each  element  of  the  array  is  processed  by  a 
quadrature  hybrid  (QH).  A  QH  is  a  power  divider  which  separates  its 
input  signal  into  quadrature  components.  We  will  refer  to  one  QH 
output  signal  as  the  in-phase  component,  Xj-j(t),  and  the  other  as 
the  quadrature  component,  Xq^(t).  The  subscript  i  indicates  the 
antenna  element  with  which  the  signal  is  associated.  We  assume  N 
antenna  elements  and  therefore  i=l,2,-,N.  The  signals  are  used 
to  define  a  signal  vector 
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Figure  1 1 - 1 .  The  IMS  array. 


The  signal  vector  is  subscripted  with  an  r  to  indicate  a  vector  of 
real  components.  Later,  when  we  introduce  analytic  signals,  the  sig¬ 
nal  vector  will  have  complex  components. 


Each  of  the  QH  output  signals  is  multiplied  by  a  weight.  Let 
X  j  f  ( t )  be  multiplied  by  weight  wji  and  Xq-j  ( t )  be  multiplied  by  weight 
wg-j .  The  weights  are  used  to  define  a  weight  vector, 


(11-2) 


vw 

The  weighted  signals  are  summed  to  produce  the  array  output,  s(t): 


s(t)  =  wrTxr(t)  =  xrT(t)wr 


( 1 1-3) 


The  IMS  array  requires  a  reference  signal,  r(t)  in  Figure  1 1  - 1 , 
which  is  an  estimate  of  the  incoming  desired  signal.  The  array  out¬ 
put  is  subtracted  from  the  reference  signal  to  form  an  error  signal, 


c(t)  =  r(t)  -  s ( t ) 


(H-4) 


The  feedback  network  adjusts  the  weights  to  minimize  the  mean  square 
value  of  the  error  signal,  E{e^(t)}  where  E{ • }  denotes  the  expecta¬ 
tion. 

The  weight  vector  that  minimizes  the  mean  square  error  satifies 
the  system  of  equations, 


vWr  E{ r  (t) 1  =  0  , 


( 1 1-5) 


where  v^„  denotes  the  gradient  with  respect  to  the  weight  vector,  Wr. 
Substituting  Equations  ( 1 1  -  3 )  and  ( 1 1 -4 )  into  Equation  ( 1 1  -  5 )  yields. 


vw  [r2(t)-2WrTE{Xr(t)r(t)l+WrTE{Xr(t)XrT(t))Wr]  =  0  .  (11-6) 
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The  equation  may  be  simplified  be  defining  the  Covariance  Matrix, 

V  =  E{Xr(t)XrT(t)}  ,  (11-7) 

and  the  Reference  Correlation  Vector, 


Sr  =  E(Xp( t)r(t) } 


( 1 1-8) 


Using  these  definitions,  equation  ( 1 1 - 6 )  becomes, 


7u  [r?(t)-2W  TS  +W  Tt  W  ]  =  0 
wrL  '  '  r  r  r  r  r 


(11-9) 


Taking  the  gradient  as  indicated  gives. 


-2S  +  2*  W  =0  .  (11-10) 

r  r  ropt 

Therefore,  the  optimum  weights  with  respect  to  minimum  mean  square 
error  are  given  by. 


W 


ropt 


(11-11) 


To  make  the  adaptive  array,  we  need  a  feedback  network  which 
will  make  the  weights  converge  to  Wr  j.-  In  this  work,  we  assume  the 
weights  are  controlled  by  the  LMS  algorithmic  ].  This  algorithm 
implements  the  gradient  technique  given  by. 


dwn 

dt 


(11-12) 


and 


dwQi 

dt 


e2(t) 


(11-13) 


where  k  is 


a  gain  constant. 


In  practice,  k  must  be  chosen  small  enough  that  Wj^  and  wg^  do  not 
respond  to  the  instantaneous  fluctuations  in  c^(t),  but  only  to  the 
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overage  value  of  r^(t)[7].  Hence  we  may  approximate  Equations  (11-1?) 
and  (11-13)  by  introducing  expected  values. 


and 


dw 


Ii 

dt 


■k  i  :EU2(t)) 
I  1 


dw, 


Qi 

dt 


■k  ^-Eu’lt)) 


(11-14) 


(11-15) 


Equations  (11-14)  and  (11-15)  can  be  combined  with  vector  notation 
into 


dW 

r 

dt 


-k  v  E{e2(t) } 


(11-16) 


Taking  the  gradient  as  was  done  for  Equation  (II-9)  gives 


dW 

dtr  =  -k[-2Sr+2VHr]  , 


(11-17) 


or 


dW 

dt 


--  +  2k  «  W 
r  r 


2k  S 


(11-18) 


In  steady  state,  when  dWr/dt=0,  we  have. 


2k  *rWr  =  2k  Sr 
ss 


(11-19) 


or 


Wr  =  Sr  •  (11-20) 

ss 

We  see  that  the  steady  state  weights  in  the  LMS  algorithm  equal  the 
optimum  weights  calculated  in  Equation  (11-14) 

The  feedback  circuit  that  implements  the  LMS  algorithm  is  easily 
found.  Substituting  Equations  ( 1 1 -4 )  into  Equations  (II-l?)  yields 
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(11-21) 


dt1  =  -k  3w-r  [r2(t)-2r(t)s(t)+s2(t)] 

wh i c h  becomes, 
dw. . 

d{1  =  -k  [-2r(t)Xn(t)-2  Xj.(t)s(t)] 

Factoring  out  the  Xj.(t)  term,  we  have, 
dw.  ■ 

d{l=  2k  Xii(t)t(t)  • 


(11-22) 


(11-23) 


A  similar  result  is  found  for  the  quadrature  weights. 


2k  XQ.(t)f(t) 


(11-24) 


From  either  Equation  (11-23)  or  Equation  (11-24),  we  find  the  feedback 
circuit  of  Figure  1 1 -2 .  A  feedback  circuit  of  this  type  is  required 
for  each  weight  in  the  adaptive  array. 


C.  Analytic  Signal  Notation 

As  was  noted  earlier,  it  is  useful  to  describe  the  LMS  array 
with  analytic  signal  notation.  To  do  this,  we  observe  that  the  two 
outputs  of  a  quadrature  hybrid  are  related  by  the  equation. 


XIi  =  XQi  ’  X0i  "XIi 
where  indicates  the  Hilbert  transform[8] , 


(H-25) 


9(0  -  f 

-  oo 


(11-26) 


Hence  we  may  combine  the  two  outputs  of  each  quadrature  hybrid 
into  an  analytic  signal. 


X^t)  = 


'Ii 


+  J  XQi  =  XTi  +  j  X 


Ii 


'Ii 


(11=27) 
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figure  1 1 -2.  The  LMS  array  feedback  loop. 


Wo  also  combine  each  pair  of  quadrature  weights  into  a  single  com¬ 
plex  weight. 


'Ii 


-  J  w 


Qi 


(11-28) 


A  minus  sign  is  used  here  to  reduce  the  number  of  conjugates  required 
in  the  resulting  equations.  We  also  define  the  analytic  signals. 


r(t)  =  r(t)  +  j  r(t) 


(11-29) 


c'(t)  =  c  ( t )  +  j  c(t) 

and 


(11-30) 


$(t)  =  s  ( t )  +  j  s  ( t ) 


(H-31) 
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for  the  reference  signal,  the  error  signal  and  the  array  output. 


The  analytic  signals  X.(t)  may  be  combined  into  a  complex  signal 
vector  X,  1 


(11-32) 


Similarly,  the  complex  weight  vector  is  defined  to  be 


(11-33) 


By  writing  out  the  real  and  imaginary  parts,  one  can  verify  that  the 
analytic  signal  corresponding  to  the  array  output  is  given  by 


s(t)  WTX  -  XTW 


(11-34) 


The  differential  equation  for  the  LMS  array  in  terms  of  analytic 
signals  is  easily  derived.  First  we  combine  Equations  (11-23)  and 
(11-24)  into  the  complex  equation 

dwn  dwn. 

dJ  -  J  dJl  2k[ElXIi(t),  (t)l  -  j  E{Xp.(t)c(t))]  .(11-35) 


Using  the  equality  [9] 


E{Xn(t)r(t)l=  E{  Xj  .(t)e(t)}  =  E{  XQj  (t)e(t) } 


(11-36) 


and  the  similar  relation  for  E{ Xn .( t)c(t) ) ,  we  can  rewrite  Eq.  (II-  35) 


dt  (WIi  -  jV  =  »cEf Xji{t)c(t)  +  Xqi(t)c(t) 


jXQ.(t)c(t)  +  jXj.(t)r(t)  } 


(H-37) 
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Now  recognizing  that 


Et  Xn (t)r(t)  +  XQ1(t)e(t)  -  jXQi(t)e(t)  +  jXJ|(t)e(t)>  = 
E{X*(t)c(t)}  ( 

we  find 
dw. 

dt-  =  k  E{X.*(t)c(t)>  .  C 

In  vector  notation.  Equation  (11-39)  becomes 

"Y  =  k  El X*(t )e(t) }  .  (] 

If  we  replace  r(t)  in  Equation  (11-40)  with 

c(t)  =  r(t)-s(t)  =  r(t)-XTW  ,  (I 

we  find  that, 

^  =  k  [E{ X*(t)r(t) }  -  E { X* ( t ) XT ( t ) }W]  .  (I 

In  this  complex  notation,  we  define  the  covariance  matrix  to  be 
*  =  El X*( t )XT( t ) }  ,  (I 

and  the  reference  correlation  vector  to  be 

S  =  El X*( t)r(t) }  .  (I 

These  definitions  allow  us  to  write  Equation  (11-42)  in  the  form, 
+  k't'W  =  kS  (I 

The  complex  steady  state  weights,  found  when  dW/dt^O  are  given  by 


r 


Wss  -  r]S  .  (11-46) 

It  is  sometimes  convenient  to  break  the  signal  vector  into 
three  components. 


X  '  <d  *  V  *n  •  (II-47) 

where  Xd  is  due  to  desired  signal,  X-j  is  due  to  interference,  and  Xn 
is  due  to  noise.  If  we  assume  that  the  three  components  are  uncorre¬ 
lated,  we  find 


<1  =  Ei  X*XT }  =  E{X  *X.T)  +  E{X,*X.T}  +  E{X  *X  T) 
d  d  it  n  n 


Using  the  notation. 


and 


"d =  EfVxd  } 


=  E 1  X . *X .  } 


*n  =  E'VY} 


we  can  write  the  covariance  matrix. 


(11-48) 


(11-49) 

(11-50) 

(H-51) 


4-  =  4>d  +  4>i  +  4>n  •  (II-52) 

We  see  that  dividing  the  signal  vector  into  three  uncorrelated  com¬ 
ponents  allows  us  to  write  the  covariance  matrix  as  the  sum  of  three 
matrices,  each  involving  one  signal  only.  This  notation  will  be  use¬ 
ful  in  Chapter  V. 

The  analytic  signal  representation  of  the  LMS  array  can  be  visual 
ized  in  terms  of  a  block  diagram.  Rewriting  Equation  (11-39)  without 
the  expectation  operator,  we  have, 

dw. 

dtl=  kX.*(t)e(t)  .  (H-53) 


This  equation  corresponds  to  the  feedback  circuit  of  Figure  1 1 - 3 . 


SIGNALS  FROM 
OTHER  ELEMENTS 


Figure  11-3.  Complex  LMS  loop. 


The  advantage  of  the  analytic  signal  notation  can  be  seen  by 
comparing  Equations  (11-20)  and  (11-46).  To  calculate  steady  state 
weights  with  real  notation,  one  must  find  the  inverse  of  iy,  a  2N 
by  2N  matrix  for  an  N  element  array.  However,  in  the  analytic  signal 
notation,  is  an  NxN  matrix  and  its  inverse  is  much  easier  to  compute. 


D .  Summa  rj/ 

The  LMS  adaptive  array  was  introduced  in  this  chapter  as  back¬ 
ground  for  the  following  chapters.  We  noted  that  the  LMS  array  mini¬ 
mizes  the  mean-square  error  between  the  array  output  and  a  so-called 
reference  signal.  In  the  remaining  chapters,  we  address  the  problem 
of  how  to  obtain  a  suitable  reference  signal  in  a  binary  FSK  communi¬ 
cation  system. 
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CHAPTER  III 

PROPERTIES  OF  THE  DESIRED  SIGNAL  BIT  STREAM 


A.  Introduction 

The  IMS  array  compares  the  array  output  signal  to  a  reference 
signal  and  minimizes  the  mean  square  error.  He  now  turn  to  the  pro¬ 
blem  of  obtaining  a  reference  signal  in  a  binary  FSK  communication 
system. 

The  technique  we  will  study  assumes  the  transmitted  bits  in 
the  FSK  signal  come  from  a  first-order  Markov  source[10].  We  will 
take  advantage  of  the  Markov  bit  transition  probabilities  to  ma^e 
a  one  bit  prediction  of  the  desired  signal  that  can  be  used  to  form 
the  reference  signal. 

In  this  chapter,  we  will  describe  the  overall  system  concept 
and  present  a  method  of  modifying  the  transmitted  bit  stream  to  con¬ 
vert  it  to  a  Markov  source.  In  Chapter  IV,  we  will  show  how  a  bit 
prediction  scheme  may  be  used  to  form  a  reference  signal  for  the 
adaptive  array. 


B.  The  Desired  Signa1_ 

Figure  1 1 1  - 1  shows  a  block  diagram  of  the  conmunicati on  system 
envisioned  in  this  study.  A  source  of  binary  information  is  assumed 
to  be  given.  The  information  bits  from  this  source  are  sent  to  a 
receiver  equipped  with  an  adaptive  array.  As  shown  in  Figure  1 1 1  - 1 , 
the  information  bits  will  first  be  coded  into  a  new  set  of  bits  and 
then  the  new  coded  bits  will  be  transmitted  to  the  adaptive  array  and 
receiver.  At  the  receiver,  the  array  output  is  first  demodulated  to 
regain  the  coded  bits,  and  these  are  then  passed  through  a  decoder 
to  recover  the  original  information  bits. 

The  purpose  of  coding  the  transmitted  bits  is  to  cause  them  to 
have  the  properties  of  a  Markov  source[10].  These  properties  allow 
the  generation  of  a  reference  signal  at  the  receiver.  The  coding 
which  converts  the  information  bits  to  a  Markov  source  will  enable  us 
to  make  bit  predictions  at  the  adaptive  array.  Bit  predictions  are 
necessary  because  Lhe  reference  signal  in  the  adaptive  array  must  be 
a  real-time  estimate  of  the  desired  signal.  Since  a  received  bit  is 
not  detected  until  the  end  of  its  bit  interval,  knowledge  of  this  bit 
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is  not  available  until  it  is  too  late  for  use  in  the  reference  signal. 
However,  if  each  bit  is  used  to  predict  the  next  bit  to  be  received 
(i.e.,  one  bit  in  advance),  an  estimate  of  the  desired  signal  can  be 
available  in  time  for  use  in  the  reference  signal.  The  purpose  of 
coding  the  information  bits  into  a  Markov  source  is  to  make  it  poss¬ 
ible  to  make  such  a  one  bit  prediction,  with  a  better  than  average 
probability  of  being  correct.  Of  course,  if  useful  information  is  to 
be  transmitted  and  received,  we  cannot  predict  every  bit  correctly. 

The  coding  process  will  allow  us  to  keep  the  number  of  bit  predic¬ 
tion  errors  low  enough  to  produce  a  suitable  reference  signal  in  the 
FSK  system. 


C.  A  Markov  Source 

To  describe  the  coding  procedure  for  the  FSK  bits,  we  need  to 
look  at  some  statistics  of  the  FSK  bit  stream.  Therefore,  we  will 
introduce  notation  for  the  bit  probabilities  :’n(0)  will  indicate 
the  probability  of  a  0  during  the  n^  Ml  '  ■  •  .al .  Likewise,  Pn(l) 

will  be  the  probability  of  a  1  on  the  -t1  -  The  notation  Pn ( 0/ 1 ) , 
will  denote  the  probability  of  a  0  or  ■"•••■  ;  nit  interval  given  that 
a  1  occured  on  the  (n-l)th  bit.  The  or-1  w  '  .ve  conditional  probabil¬ 
ities,  Pn(0/0),  Pn(l/1)  and  Pn(l/0),  a>  ■■  •'  :nen  in  the  same  manner. 

The  coding  procedure  that  we  will  use  at  the  transmitter  con¬ 
verts  the  information  bit  sequence  into  a  Markov  bit  sequence[l 0] . 

In  an  nth  order  Markov  sequence,  the  ivobability  of  the  next  bit  is 
a  function  of  the  previous  n  bits.  We  will  consider  in  particular, 
a  first  order  Markov  bit  sequence.  For  this  case,  the  probability 
of  a  1  or  0  on  each  bit  interval  depends  only  on  which  bit  was  trans¬ 
mitted  previously.  This  relationship  may  be  visualized  in  terms  of 
the  transition  diagram  as  shown  in  Figure  1 1 1-2.  The  arrows  in  the 
diagram  represent  bit  transitions  from  the  ( n - 1 ) t n  5^  interval  to 
the  n'-'1  bit  interval.  We  can  add  the  probabilities  and  write  the 
equations. 


P„(0)  =  Pn(0/0)Pn_1(0)  +  Pn(0/l)Pn--|(l)  (iii-D 

and 

Pn(D  =  Pn(i/o)Pn1(o)  +  PR (1/1 )Pn_ i ( 1 )  •  (111’2) 

We  will  further  assume  that  the  first  order  Markov  sequence  is 
stationary  and  symmetric.  By  stationary,  we  mean  that  the  probabi¬ 
lities  do  not  change  from  bit  interval  to  bit  interval.  Therefore, 
we  can  drop  the  subscripts  from  Equations  ( 1 1 1  - 1 )  and  (III-2)  and 
wr  i te , 
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0 

BIT 


Pn  (  0/0) 


0 

BIT 


!"  kju  re  1 1 1  -  2 .  Transition  diagram  for  a  first-order 
Markov  source. 


P(0)  -  P(0/0)P(0)  +  P(0/1)P(1)  (111-3) 

and 

P(l)  3  P(1/0)P(0)  +  P( 1/1 )P( 1 )  .  ( 1 1 1 -4) 

The  assumption  that  the  Markov  sequence  is  symmetric  means, 

P(0/1)  =  P(l/0)  =  q  (HI-5) 

and 

P(0/0)  -  P(l/1)  =  1-q  .  (111-6) 

Substituting  Equations(llI -5)  and  ( 1 1 1 -6)  into  Equations  (111-3) 
and  (111-4)  gives. 
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P(0)  ~  P(l) 


1 

2 


( 1 1 1-7) 


In  summary,  a  first  order  stationary  and  symmetric  Markov  bit 
sequence  is  described  by  the  equations, 

P(0)  =  P(l)  =  1  ,  (HI-8) 

P(bit  change)  =  q  ,  ( 1 1 1  - 9 ) 

and 

P(no  bit  change)  =  1-q  .  (III-10) 


At  the  receiver,  the  Markov  nature  of  the  desired  signal  allows 
us  to  make  bit  predictions.  Assume  for  example  that  the  probability 
of  a  bit  change  is  greater  than  one  half. 


P(l/0)  =  P(0/1)  >  - 


(111-11) 


Then  we  can  predict  the  next  bit  by  assuming  it  to  be  the  opposite  of 
the  last  bit.  The  probability  of  a  correct  prediction  is  equal  to  q, 
the  bit  change  probability.  The  probability  of  a  bit  prediction  error 
is  equal  to  1-q. 


D.  Coding  and  Decod incj  of  the 

Message  Bits 

As  noted  above,  we  want  to  transmit  a  Markov  bit  stream  to  the 
adaptive  array.  Therefore  to  send  an  arbitrary  binary  message,  we 
must  first  convert  the  message  bits  into  Markov  bits.  We  will  now 
illustrate,  in  the  form  of  an  example,  how  an  arbitrary  bit  sequence 
may  be  converted  into  a  Markov  bit  sequence. 

A  method  of  forming  a  Markov  bit  stream  from  an  arbitrary  bit 
sequence  with  equally  likely  bit  transitions  is  illustrated  in  Figure 
II 1-3.  In  this  illustration,  each  group  of  3  message  bits  is  mapped 
into  a  group  of  four  Markov  bits  for  transmission  to  the  adaptive 
array.  The  mapping  is  done  by  first  listing  all  possible  4  bit  words 
and  counting  the  number  of  bit  transitions  in  each  word.  Then  we 
choose  the  8  words  which  have  the  highest  number  of  bit  transitions. 
The  set  of  8  words  chosen  are  marked  with  a  check  in  the  figure.  In 
figure  III -3( b)  we  list  all  possible  combinations  of  3  bits  which  are 
the  message  words.  Then  to  each  of  these  3  bit  message  words,  we 
assign  one  of  the  4  bit  words  marked  with  a  check.  These  assignments 
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4  8 IT 
WORDS 


NUMBER  OF 
TRANSITIONS 


OOOO 

0 

0001 

1 

MESSAGE 

TRANSMITTED 

0010 

2  / 

BITS 

BITS 

00  1  1 

1 

000 

00  1  0 

0  100 

2  / 

001 

0  100 

0  101 

3  y 

| . 

0  101 

0  i  10 

2  V 

01  1 

0  110 

0  1  II 

1 

100 

1  001 

1000 

1 

101 

10  10 

1001 

2  y 

1  10 

1011 

1010 

3  y 

1  1 1 

110  1 

1011 

2  y 

(b) 

1  1  00 

i 

1  101 

2  y 

1  1  10 

i 

1  1  1  l 

0 

(a) 


Figure  1 1 1 - 3 .  Markov  coding  example. 


19 


•ire  used  at  the  transmi t  ter  so  that  for  each  3  message  bits,  4  Markov 
hits  are  transmitted.  We  find  for  the  Markov  bits. 


<1  -  P { bit  transition!  -  0.69 


(311-12) 


We  note  that  by  transmitting  the  Markov  bits  instead  of  the 
message  bits,  the  number  of  transmitted  bits  is  increased  by  the 
ratio  4/3.  If  the  message  bit  rate  is  held  constant,  the  coding  pro¬ 
cess  will  require  an  increased  transmission  bit  rate.  Consequently, 
the  bandwidth  of  the  FSK  siynal  must  be  increased  as  well. 

The  bandwidth  increase  can  he  related  to  the  transition  proba¬ 
bilities  of  the  Markov  bits.  Let  denote  the  entropy  of  the  message 

bits  and  Hq  the  entropy  of  the  coded  Markov  bits.  The  coded  bit 
stream  is  a  symmetric  first  order  Markov  sequence  and  therefore[ 1 1 ] 

HC  =  9  1o92  J  +  0-9)  log2  ■  (111-13) 


let  R^  denote  the  bandwidth  of  the  FSK  signal  when  the  message  bits 
are  transmitted  without  being  coded  into  Markov  bits.  Bq  will  be  the 
bandwidth  of  the  FSK  signal  when  the  Markov  bits  are  transmitted. 

Since  the  message  bit  rate  is  held  constant,  we  can  write  the  required 
bandwidth  of  the  transmitted  Markov  bits  as 


Br 


--  R, 


H 

m  h; 


M 


(311-14) 


'min  C 

Also  we  can  define  the  minimum  bandwidth  spreading  factor, 
H. 


F  = 


M 

h; 


(111-15) 


If  we  assume  that  the  message  bits  have  unit  entropy. 


H 


M 


1 


we  can  write. 


( 1 11-16) 


F 


1 

q  1°92  J  +  (!-q)  log2 


(111-17) 
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A  ■->  1  o  t  of  Equation  (III-17)  is  shown  in  Figure  1 1 1  -4 .  We  note  that 
as  a  result  of  doubling  the  bandwidth  of  the  transmitted  FSK  signal, 
we  can  make  bit  predictions  with  a  0.9  probability  of  being  correct. 

Equation  (III-17)  also  allows  us  to  compare  the  coding  example 
of  Figure  1II-3  to  theoretical  results.  In  the  example,  a  bandwidth 
spreading  ratio  of  4/3  resulted  in  a  bit  change  probability  of  0.7. 
From  the  plot  in  Figure  III -4 ,  we  see  that  the  maximum  theoretical 
bit  change  probability  for  the  example  case  is  0.8.  Although  the 
coding  scheme  in  the  example  does  not  give  optimum  results,  it  does 
demonstrate  the  feasibility  of  the  idea.  In  this  report,  we  will 
not  consider  the  general  problem  of  converting  a  bit  stream  into  a 
Markov  bit  stream  of  given  characteristics.  Rather,  we  will  just 
use  this  example  to  illustrate  that  it  can  be  done. 


0.5  0.6  0.7  0.8  0.9  1.0  q 


Figure  III -4 .  Bandwidth  spread  versus  prediction 
probabil i ty. 
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CHAPTER  IV 

REFERENCE  SIGNAL  GENERATION 


A.  Introduction 

The  previous  chapter  presented  a  coding  procedure  for  the  bit 
stream.  The  bits  are  coded  to  make  the  bit  stream  first-order  Markov, 
by  detecting  the  last  bit  received  and  knowing  the  bit  transition 
probabilities  of  the  Markov  code,  we  can  make  bit  predictions  at  the 
adaptive  array.  In  this  chapter,  we  will  describe  the  operation  of 
the  adaptive  array  reference  signal  generation  circuity  and  will  show 
how  the  bit  predictions  are  used  in  forming  a  reference  signal  from 
the  array  output.  Next  we  will  demonstrate  the  formation  of  a  refer¬ 
ence  signal  when  the  adaptive  array  receives  CW  interference  in 
addition  to  desired  signal.  The  examples  in  this  chapter  will  assume 
that  perfect  bit  predictions  are  made  at  the  receiver.  In  the  next 
chapter,  we  will  look  at  the  effects  of  bit  prediction  errors  on  the 
performance  of  the  system. 


B.  The  Adaptive  Array  and  the_ 

FSK  Detector 

Figure  IV - 1  shows  a  block  diagram  of  an  adaptive  array  with  its 
reference  signal  generation  circuitry.  The  array  is  the  LMS  array 
described  in  Chapter  II.  An  incoherent  FSK  bit  detector  is  connected 
to  the  array  output  and  produces  the  received  bit  stream.  A  pre¬ 
diction  circuit  looks  at  the  received  bits  and  makes  bit  predictions 
according  to  the  Markov  properties  of  the  bits.  The  array  output  and 
the  bit  predictions  are  used  by  the  reference  signal  processing  loop 
to  produce  a  reference  signal. 

The  detector  behind  the  LMS  adaptive  array  is  assumed  to  be 
incoherent.  The  block  diagram  of  an  incoherent  FSK  detector[12]  is 
shown  in  Figure  IV-2.  The  signal  is  split  two  channels,  each  contain¬ 
ing  a  bandpass  filter  and  an  envelope  detector.  In  the  upper  or  0 
channel,  the  filter  has  center  frequency  u0  corresponding  to  the  0 
frequency  of  the  FSK  signal.  The  filter  in  the  lower  channel  has 
center  frequency  ^i  corresponding  to  the  1  frequency.  The  output 
of  each  filter  is  then  envelope  detected.  If,  for  example,  at  the 
end  of  each  bit  interval,  the  output  of  the  envelope  detector  in  the 
0  channel  is  greater  in  magnitude  than  that  in  the  1  channel,  a  0  bit 
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ARRAY 

OUTPUT 


Figure  I V - 2 .  Incoherent  FSK  bit  detector. 


selected.  Likewise,  when  the  output  of  the  1  channel  is  greater  a 
1  iiit  i  s  detec  ted  . 

During  a  particular  bit  interval,  one  channel  of  the  FSK  detec¬ 
tor  contains  both  signal  and  noise  while  the  other  channel  contains 
only  noise.  Bit  detection  errors  are  made  whenever  the  magnitude  of 
the  noise  in  the  noise  only  channel  exceeds  the  magnitude  of  the 
signal  plus  noise  in  the  other  channel.  We  will  assume  temporarily 
in  this  section  that  no  bit  detection  errors  are  made.  Later,  we 
will  look  at  the  effects  of  bit  detection  errors  on  the  performance 
of  the  system. 


C.  The  Reference  Signal  Processing  Loop 

A  block  diagram  of  the  reference  signal  processing  loop  is  shown 
in  Figure  I V -3 .  The  array  output  and  the  bit  predictions  are  inputs 
to  the  loop  and  the  loop  output  is  the  reference  signal.  In  the  fol¬ 
lowing  discussion,  we  will  look  individually  at  the  bandpass  filter, 
the  local  oscillator,  and  the  limiter.  Then  we  will  show  how  the  loop 
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MIXER  1 


operates  when  the  array  output  consists  of  both  desired  signal  and 
interference. 

The  bandpass  filter  is  used  in  the  loop  to  separate  the  desired 
and  the  undesired  components  of  the  array  output.  A  typical  ampli¬ 
tude  response  of  the  filter  is  illstrated  in  Figure  IV -4 .  The  filter 
has  a  center  frequency  of  u>c  and  its  half  power  bandwidth  is  B. 


Figure  IV -4 .  Transfer  function  of  the  reference 
loop  bandpass  filter. 


The  oscillator  shifts  the  frequency  of  the  array  output  to  the 
center  of  the  bandpass  filter.  Ideally,  if  wosc  is  the  frequency  of 
the  oscillator  signal, 

u>  -u>  during  a  0  bit 

"'osc  =  4  (IV-1) 

u>c-<.ii  during  a  1  bit 

We  do  not  know  the  bit  being  received,  however,  until  the  end  of  the 
bit  interval.  Therefore,  we  set  the  oscillator  frequency  according 
to  the  predicted  bit  for  the  current  bit  interval.  The  oscillator 
frequency  then  is 
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{uic-(oq  when  a  0  is  predicted 

w  -w-j  when  a  1  is  predicted 


( IV- 2  J 


A  limiter  is  included  in  the  loop  so  the  reference  signal  will 
have  fixed  amplitude.  A  typical  transfer  function  for  a  limiter  is 
shown  in  Figure  I V - 5 .  Below  the  saturation  voltage  R,  the  device  be¬ 
haves  like  a  linear  amplifier  with  a  gain  of  G.  When  the  device  is 
saturated,  the  output  signal  is  clipped  and  remains  below  voltage  R. 


Figure  IV- 5 .  Transfer  function  of  the  limiter. 


D.  Operation  of  the  Reference 
Signal  Loop 

We  will  assume  in  our  discussion  of  the  reference  signal  loop 
that  the  array  receives  both  an  FSK  desired  signal  and  a  CW  interfer¬ 
ence  signal.  The  incoming  desired  signal  at  the  array  elements,  in 
real  notation,  is  given  by, 


d  (t)  =  /2  Ad  cos(o>dt+*d)  ,  ( I V - 3 ) 
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whore 


during  a  0  bit 
during  a  1  bit 


(IV -4 ) 


We  use  /2  Ad  for  the  signal  amplitude  so  that  the  outputs  of  the  quad¬ 
rature  hybrids  will  have  amplitude  Ad.  The  phase  angle,  vd »  is  assum¬ 
ed  uniformly  distributed, 

r  1 

A  -  -1I<W  .<n 

2’  d  .  ( IV- 5 ) 

0  elsewhere 


pUd) 


The  desired  signal  component  of  the  array  output  will  be  of  the  form. 


sdU)  =  Bd  cos(u)dt+Yd)  ,  ( I V - 6 ) 

where  the  amplitude,  B,,  and  the  phase,  yd, incorporate  the  gain  and 
phase  shift  caused  by  the  array  weights.  Similarly,  we  assume  the 
interference  signal  is  given  by, 

i  (t)  =  /2  Ai  cos(a-it+^.)  .  ( IV -7 ) 


The  magnitude  of  the  interference  is  /2  A.  and  u>.  is  the  interference 
frequency.  The  phase  angle,  ,  is  assumed  uniformly  distributed. 


p(i'i)  = 


-n<4)^<n 

el sewhere 


(IV -8 ) 


The  array  output  component  associated  with  the  interference  is. 


Sj(t)  =  COs(u)ft+y.) 


( IV  -  9 ) 


As  with  the  desired  signal,  the  amplitude  B-j  and  the  phase  angle  Y-j 
indicate  the  effects  of  the  array  weights.  The  total  array  output 
due  to  desired  signal,  interference,  and  noise  is, 
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s(t)  =  sd(t)+si(t)+sn(t) 


(IV-10) 


We  will  assume  that  the  desired  signal  to  noise  ratio  is  high  enough 
to  neglect  noise  in  the  reference  loop.  This  assumption,  however, 
does  not  imply  that  noise  is  not  a  factor  in  the  bit  detection  process 
in  the  receiver.  Neglecting  noise,  we  can  write  the  array  output. 


S  ( t )  =  Bd  COs(a>dt+Yd)  +Bi  COS(u.it+Yi  )  -  (IV-11) 

Refering  to  Figure  I V - 3 ,  we  see  that  the  array  output  enters 
the  reference  signal  loop  and  is  mixed  with  the  oscillator  signal. 

For  now,  we  will  assume  that  perfect  bit  predictions  are  made  and 
therefore  the  oscillator  frequency  is 


“'osc  =  v^d  •  u 

The  problems  caused  by  bit  prediction  errors  will  be  discussed  in 
Chapter  V. 

The  output  of  the  loop  oscillator  is  given  by, 


y(t)  =  cos[(V(od)t+Yosc] 


( I V  - 1 3 ) 


where  the  amplitude  is  assumed  to  be  unity  and  Yosc  the  phase  angle. 

The  output  of  Mixer  1  is  the  product  of  the  oscillator  signal 
and  the  array  output.  The  mixer  output  is 


f(t)  =  [Bd  cos(W(jt+Yd)+Bi  cos(u).t+Yi)]  C0S[(v‘Vt+>0SC] 

( I V  - 1 4  ) 

equation  (IV-14)  can  be  expanded  into 

f(t)  --  /  COs[mct+Yd+Y()SC]  +  2d  COs[((l,c-2,od)t-Yd+Vosc] 

B .  B . 

+  0--  COs[(‘n  +m.  -(0  ,  )t+Y  -+Y  ]  +  o'1"  COs[(ui  -<■! . -m  . )  t- y  ■  f  r  ]■ 

2  ci  dy  i  osc  2  L '  c  i  d  i  'osc 

( IV- 15) 
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The  amplitude  response,  H(u>),  and  the  phase  response,  e(,.),  of 
the  bandpass  filter  are  shown  in  Figure  I V -4 .  We  will  assume  that  the 
frequencies  uc-9  and  wc- (m-j+u'd )  are  outside  the  bandpass  of  the 
filter.  Therefore,  we  can  drop  these  components  from  the  filter  input 
and  write. 


L  B. 

f  *  (t)  r  2  -  CI,sUct^V \scj  4  2  c  0  5  f.  (  wc  J'W  .•  -  )  t+\  i  +  >0  sc  ] 


( IV- 16) 


The  filter  output  signal  is  given  by, 

g(t)  =  2J-U{,.,c)cos(..cttvd+w«(„c)) 
B. 


+  H(io  +io  . -to  ,  )cOS  (  (to  +(J  . -w  .  )  t+y  .+y  +  o(u>  +u.-io. )) 

2  vc  l  d'  "  c  i  d  '  i  osc  c  l  d 

( 1 V  - 1  7 ) 

At  Mixer  2,  the  output  of  the  bandpass  filter  is  multiplied  by 
the  oscillator  signal.  The  second  mixer  output  is  given  by, 

p(t)  =  j/  H((„c)cOs[mct+Td+yosc  +  0(mc)] 

Bi  1 

4  2' '  H  ( lllc  +("  i  ‘  u’d  ^ c  0  s  ^  wc  +l°  i  ‘ wd  ^ 1 +T  i +  Yo  sc + 0  ^c+li'  i  ”  ‘"’d  ^  [ 

COS[(()c+a'd  )t+Y0SC]  ■  (  I V  - 1 8 ) 


We  will  neglect  the  terms  in  Equation  ( IV- 1 8 )  with  frequencies  of  the 
order  of  2<0(j  under  the  assumption  that  the  transient  response  of  the 
adaptive  array  is  too  slow  to  respond  to  them.  If  this  assumption 
were  not  valid,  we  would  add  a  zonal  filter  at  the  output  of  Mixer  2. 
In  either  case,  we  can  drop  the  high  frequency  terms  and  write  the 
second  mixer  output. 


p(t)  --  4d  H(u.c)cos[a>dt  +  Yd+O(u>c)] 


H(to  +io  . -dt  ,  )cOs[io  •  t+y .+0  (to  +ut . -to  , )  ] 


(IV-19) 


The  adaptive  array,  in  the  absence  of  noise,  will  make  the  mag¬ 
nitude  of  the  array  output  track  the  magnitude  of  the  reference  signal. 
If  the  reference  signal  amplitude  is  greater  than  the  array  ojto"t 
amplitude,  the  weights  and  the  array  output  will  continuously  .  rease. 
A  limiter  is  included  in  the  reference  signal  loop  to  keep  the  array 
output  at  a  fixed  level,  inside  the  dynamic  range  of  the  electronics. 
The  limiter  is  connected  to  the  output  of  Mixer  2,  and  its  input  is 
given  by  Equation  (IV-19).  As  illustrated  in  Figure  I V - 5 ,  the  limiter 
behaves  like  a  linear  amplifier  below  saturation. 

We  will  assume  that  both  terms  in  Equation  (IV-19)  fall  below 
the  limiter  cutoff  voltage.  This  assumption  is  valid  because  of  the 
relationship  between  the  magnitudes  of  the  reference  signal  and  the 
array  output.  If  the  limiter  were  in  the  clipping  region  and  if  the 
reference  signal  magnitude  were  smaller  than  the  array  output  magni¬ 
tude,  the  array  output  would  be  decreased.  The  array  output  would 
continue  to  be  reduced  until  the  limiter  input  signal  reaches  the  knee 
of  the  transfer  function.  Therefore,  we  can  write  the  limiter  output, 
which  is  the  reference  signal,  as, 

r(t)  =  ^  G  H(u)c)cos[u»ct+Yd+e(loc)] 

B  . 

+  G  H(w^+m,. -a)^  )cos[<o^  t+y  .  +6>  (w  +ui,. )  ]  •  (  IV  —  20) 

The  reference  signal,  given  by  Equation  (IV- 20 ) ,  consists  of 
both  a  desired  and  an  undesired  component.  The  loop  gain  for  the 
desired  component,  measured  between  the  array  output  and  the  reference 
signal ,  is. 


L 


d 


(IV-21  ) 


ps 

* 


The  loop  gain  for  the  interference  component  of  the  reference  signal 
is  given  by, 

l.  =  G- H (*.+,„. -u<j)  •  ( IV - 22 ) 

Since  we  are  temporarily  considering  the  case  of  perfect  bit  pre¬ 
dictions,  the  loop  gain  for  the  desired  signal  is  not  a  function  of 
which  bit  was  transmitted.  The  loop  gain  for  the  interference  com¬ 
ponent,  however,  contains  an  uy  term  and  is  therefore  a  function  of 
the  transmitted  bit  stream.  The  average  loop  gain  for  the  interfer¬ 
ence  can  be  written. 
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tmumm 


L  .  =  P (0)  7-  H(u.  +ui ■  • 
^  I  4  c  ^ 


*0)  +  P(l)  |  H(u>c+ui-u.1) 


(IV-23) 


Tor  the  special  case  where  the  bits  are  equally  probable. 


Li  '  8  [H(u'C+u,i'°'0)  +  H(u'c+ui'u'l)] 


(IV- 24  ) 


The  relationsi.ip  between  the  amplitude  of  the  reference  signal 
and  the  amplitude  of  the  array  output  signal  is  an  important  para¬ 
meter  of  the  reference  loop.  This  parameter  is  the  loop  gain  measured 
from  the  array  output  around  the  reference  loop  to  the  reference  signal 
port.  To  explain  the  importance  of  the  reference  loop  gain,  we  will 
consider  the  case  where  the  signals  passing  through  the  loop  are  below 
the  limiter  cutoff  amplitude.  In  this  situation,  we  want  the  inter¬ 
ference  signal  passing  through  the  reference  loop  to  be  attenuated 
so  we  make  the  loop  gain  for  the  interference  less  than  one, 


L .j  "  g  [H(mc+M.-u>0)  +  H(wc+uw-io^  )] 


(IV-25) 


At  the  same  time,  we  want  desired  signal  passing  through  the  reference 
loop  to  be  amplified  so  we  require 


Ld  =  “  H(„c)  -  I 


(IV-?6) 


If  fquations  (IV-25)  and  (IV- 26)  are  satisfied,  the  desired  part  of 
the  array  output  will  increase  up  to  the  limiter  cutoff  voltage  while 
the  interference  magnitude  is  decreased.  Therefore  the  array  will 
steer  a  beam  maximum  toward  the  desired  signal  and  form  a  null  on  the 
i  nterference. 


E.  Computer  Simulation  of  the 

FSK  System 

A  computer  simulation  program  was  written  to  study  the  perform¬ 
ance  of  the  reference  signal  processing  loop.  The  program  simulates 
a  two  element  IMS  array,  a  bit  prediction  circuit,  and  a  reference 
signal  loop.  In  this  section,  we  will  describe  the  simulation  and 
present  some  results.  A  description  of  the  program  and  a  copy  of  the 
Fortran  code  are  included  in  the  Appendix  at  the  end  of  this  thesis. 

The  geometry  of  the  two  element  array  is  shown  in  Figure  I V - 6 . 
The  antennas  are  assumed  to  have  omnidirectional  element  patterns  and 
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Figure  IV-6.  The  geometry  of  the  two-element 
adaptive  array. 


to  be  separated  by  a  distance  of  L  meters.  Both  incoming  signals  are 
assumed  to  be  in  a  plane  containing  the  array  elements.  The  incident 
FSK  signal,  in  complex  notation,  is  given  by 


d(t)  =  /2  Ad 


e 


The  desired  signal  arrives  at  an  angle  measured  from  broadside. 
Since  the  frequency  of  the  desired  signal  changes  from  bit  to  bit, 
the  interelement  phase  shift  depends  on  which  bit  is  being  received. 
For  a  0  bit,  the  interelement  phase  shift  is 


♦d<-o> 


si  no 


(1 


and  for  a  1  bit. 


where  L  is  the  antenna  elanent  separation  and  c  is  the  speed  of  light 
We  assume  that  the  two  FSK  frequencies  are  close  enough  together  so 
that  we  can  average  the  two  phase  shifts  and  use' 

*d(“o)+*dh) 

-•d  ~  "2 -  UV-30: 


for  the  calculations  in  this  chapter.  The  complex  form  of  the  inter¬ 
ference  signal  notation,  is  given  by, 

j  («.  •  t+ir, ) 

1  (t)  -  /2  fli  e  1  ,  ( I V - 3 1 

and  the  interel ement  phase  shift  is  written, 


2n  L  ■  _ 

<}> .  =  -  sine  . 

N  1  w .  1 

1 


(IV-32! 


We  also  assume  that  each  received  signal  component  contains  a  noise 
term  with  power 


As  defined  in  Chapter  II,  the  signal  vector  due  to  the  desired 
signal  is 

1 

( I V  -  3  3 ' 


Xd  =  Ad 


■Jt 


The  interference  component  of  the  signal  vector  is  given  by, 


i  \ 


J  ( <X>  ^  t + ip  ) 


( IV-34 ; 


and  the  noise  component  is 


Ei  ri j  ( t )n1  * ( t ) }  =  E{n2(t)n2*(t)}  =  02  ,  ( I V - 36 ) 


and 


E{n1(t)n2*(t)}  =  0  .  (IV- 37 ) 

The  results  of  the  simulation  program  consist  of  plots  of  each 
of  the  four  weights  versus  time  and  an  antenna  pattern  taken  at  the 
end  of  the  simulation  period.  We  will  present  the  results  of  two 
simulations  in  this  section.  First  will  be  for  the  case  of  an  FSK 
desired  signal,  perfect  bit  predictions,  and  no  interference.  In  the 
second  simulation,  we  will  include  CW  interference  directly  on  one 
of  the  two  FSK  frequencies. 

A  list  of  input  data  for  the  first  simulation  example  is  shown 
in  Figure  IV -7 .  The  FSK  desired  signal  is  incident  from  broadside 
and  has  unit  magnitude.  The  selection  of  frequencies  in  the  simula¬ 
tion  is  rather  arbitrary  since  the  simulation  is  not  done  in  real  time 
Scaling  a  frequency  up  or  down  merely  changes  the  scale  on  the  time 
axis  in  the  results.  Desired  signal  frequencies  of  990  hertz  for  a  C 
and  1010  hertz  for  a  1  were  chosen  for  the  simulation.  The  center 
frequency  of  the  bandpass  filter  in  the  reference  signal  loop  was 
chosen  to  be  2000  hertz.  Also,  at  the  center  of  the  bandpass  filter, 


H(u>c)  =  1  _  (IV-  38 ) 

The  limited  gain  is  set  to  2, 

G  =  2  ,  ( I V  -  39 ) 

so  that  the  loop  gain  for  the  desired  portion  of  the  reference  signal 
is 

Ld  =  2  .  ( IV -4  0 ) 

The  antennas  are  separated  by  one-half  of  the  average  wavelength  of 
the  desired  signal.  Finally,  the  array  convergence  factor,  k,  was 
adjusted  so  that  the  weights  would  '•each  steady-state  in  approxi- 
r lately  one-half  of  the  simulation  period. 

The  results  of  the  first  simulation  example  are  presented  in 
figure  IV-8.  The  four  weight  curves  are  identical  because  the  desired 
signal  is  at  broadside  and  therefore  both  antenna  elements  receive 
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ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 


INCOMING  SIGNALS: 

DESIRED  FSK  SIGNAL 

"O"  BIT  FREQUENCY=990  HERTZ 
"1“  BIT  FREQUENCY-1010  HERTZ 
POUER -0  DB 

ANGLE  =0  DEGREES  FROM  BROADSIDE 
FIT  RAT E  =  1 0  BITS/SECOND 
MARKOV  TRANSITION  PROBABILITY^! . 

INTERFERENCE  SIGNAL 
NONE 

NOISE 

POUER /ELEMENT1- 10  DB 

ADAPTIVE  ARRAY 

ELEMENTS^ 

SEPARATIONS. 5  DESIRED  SIGNAL  UAVELENGTHS 
CONVERGENCE  F A C T OR ( K ) =0 . 00005 

REFERENCE  SIGNAL  LOOP 
OSCILLATOR 

FREQUENCY  FOR  “0"  BIT  PREDICT  1 0N= 1 01 0  HERTZ 
FREQUENCY  FOR  "1"  BIT  PREDICT I0N=??0  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCY=2000  HERTZ 
BANDUIDTHSO  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOW  CUTOFFS  DB 


Figure  IV-7.  Simulation  input  data  for  the  case  of 
no  interference. 
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the  \aine  signal.  The  steady-state  antenna  pattern  has  its  main  beam 
in  tiie  direction  of  the  desired  signal  indicating  proper  operation 
of  the  adaptive  array. 

For  the  second  simulation  example,  we  will  add  CW  interference 
directly  on  one  of  the  two  FSK  bits.  A  list  of  input  data  for  the 
second  simulation  is  shown  in  Figure  IV- 9 .  Except  for  the  addition 
of  the  interference,  all  the  input  parameters  are  unchanged  from  the 
previous  case. 

The  addition  of  the  interfernece  affects  not  only  the  adaptive 
array  response,  but  also  the  bit  detection  and  bit  prediction  pro¬ 
cesses.  During  lock-up,  the  array  output  component  due  to  interfer¬ 
ence  may  have  much  higher  power  than  the  desired  component.  As  a 
result,  the  FSK  detector  will  not  be  able  to  defect  the  desired 
signal  bits.  We  assume  in  the  simulation  program  however  that  bit 
detection  and  timing  acquisition  continue  even  when  the  desired  sig¬ 
nal  power  is  below  the  interference  power  at  the  array  output.  In 
Chapter  V,  we  will  look  at  the  problems  of  timing  acquisition  and 
bit  detection  during  lock-up. 

The  results  of  the  system  simulation  with  interference  included 
are  shown  in  Figure  I V- 1 0 .  Note  that  with  the  addition  of  interfer¬ 
ence,  the  weights  are  not  smooth  curves.  The  roughness  in  the  weight 
curves  appears  because  the  desired  signal  and  the  interference  some¬ 
times  look  alike  to  the  adaptive  array.  When  a  0  bit  is  received 
along  with  CW  interference  on  the  0  bit  frequency,  the  array  tries 
to  acquire  both  signals.  During  the  early  part  of  the  adaptation 
period,  the  tendency  to  acquire  both  signals  results  in  abrupt  changes 
in  the  weights.  As  the  null  forms  and  the  array  output  decreases, 
the  weight  curves  become  smoother.  In  steady-state,  the  interfer¬ 
ence  is  practically  eliminated  from  the  array  output  and  therefore 
the  weight  curves  become  smooth. 

The  antenna  pattern  taken  at  the  end  of  the  simulation  period 
is  also  shown  in  Figure  IV - 1 0 .  A  null  of  about  40  dB  in  the  direc¬ 
tion  of  the  interference  shows  that  the  system  is  functioning  properly. 
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ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 

INCOMING  SIGNALS: 

DESIRED  PSK  SIGNAL 

" 0 "  BIT  FREQUENCY -990  HERTZ 
"1“  BIT  FREQUENCY^  010  HERTZ 
POWERS  DB 

ANGLES  DEGREES  FROM  BROADSIDE 
BIT  F  A  T  E  =  1 0  BITS/SECOND 
MARKOV  TRANSITION  PROBABILITY^ . 

INTERFERENCE  SIGNAL 
FREQUENCYS90  HERTZ 
F'OUER-  1  0  DB 

ANGLEM5  DEGREES  FROM  BROADSIDE 
NOISE 

POUER/EIEMENT=-10  DB 

ADAPTIVE  ARRAY 

ELEMENT  S  =  2 

SEPARAT I0NS.5  DESIRED  SIGNAL  UAVELENGTHS 
CONVERGENCE  FACTOR  <  K ) =0 . 00005 

REFERENCE  SIGNAL  LOOP 
OSCILLATOR 

FREQUENCY  FOR  "0"  BIT  PR E D I C T I ON = 1 0 1 0  HERTZ 
FREQUENCY  FOR  "1”  BIT  PREDICT I0NS90  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCY=2000  HERTZ 
BANDUIDTH*1 0  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOW  CUTOFFS  DB 


icjure  IV-9.  Simulation  input  data  for  the  case 
of  CW  interference. 


CHAPTER  V 

REFERENCE  LOOP  DESIGN  FACTORS 


A.  Introduction 

In  this  chapter,  we  will  Iook  at  some  aspects  of  the  reference 
signal  loop  which  can  affect  the  performance  of  the  adaptive  array. 
First  we  will  consider  the  problem  of  bit  prediction  errors.  The 
simulation  program,  introduced  in  the  previous  chapter,  will  be  used 
to  show  the  effects  of  bit  prediction  errors  on  the  array  weights. 

The  second  problem  that  we  will  consider  occurs  when  there  is  a  phase 
shift  between  the  array  output  and  the  reference  signal.  We  will 
show  that  this  phase  error  causes  the  array  weights  to  oscillate. 
Also,  we  will  show  that  in  the  presence  of  interference,  the  signal 
to  interference  plus  noise  ratio  (SINR)  at  the  array  output  decreases 
with  increasing  phase  error.  Finally,  we  will  look  briefly  at  the 
problem  of  synchronizing  the  bit  detector  and  the  bit  predictions 
with  the  received  bit  stream. 


B.  Bit  Prediction  Errors 

We  demonstrated,  in  the  previous  chapter,  the  performance  of  the 
reference  signal  loop  for  the  case  of  perfect  bit  predictions.  If  use 
ful  information  is  to  be  received  however,  prediction  errors  cannot  be 
avoided.  Therefore,  we  are  interested  in  the  effect  of  bit  prediction 
errors  on  the  performance  of  the  FSK  system. 

A  direct  calculation  of  the  effects  of  bit  prediction  errors  on 
the  system  performance  would  be  a  difficult  task.  All  the  different 
parameters  involved  in  bit  detection  and  prediction,  reference  loop 
design  and  adaptive  array  combine  together  to  form  a  very  complicated 
problem.  Therefore  we  will  use  computer  simulation  results  to  study 
the  effects  of  bit  prediction  errors.  The  system  simulation  program 
introduced  in  the  previous  chapter  will  be  used  to  demonstrate  the 
effects  of  several  bit  prediction  error  probabilities. 

To  see  the  results  of  a  single  prediction  error,  we  modified 
the  system  simulation  program  so  a  prediction  error  would  occur  on  the 
lA^'1  bit.  A  list  of  the  program  input  data  is  given  in  Figure  V-l. 
Except  for  the  prediction  error,  the  input  data  is  the  same  as  for 
the  last  example  in  Chapter  IV.  The  FSK  desired  signal  is  incident 


ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 


INCOMING  SIGNALS: 

DESIRED  FSK  SIGNAL 

"O"  BIT  FREQUENCY=?90  HERTZ 
“I"  BIT  FREQUENCY11 1  01  0  HERTZ 
POWERS  DB 

ANGLE  =  0  DEGREES  FROM  BROADSIDE 
BIT  RAT E  =  1 0  BITS/SECOND 
MARKOV  TRANSITION  F'RQBAB IL I TY  =  1  . 

INTERFERENCE  SIGNAL 
NONE 

NOISE 

P0UER/ELENENT=-1 0  DB 

ADAPTIVE  ARRAY 

ELENENTS=2 

SEPARATIONS. 5  DESIRED  SIGNAL  UAVELENGTHS 
CONVERGENCE  F ACTOR (K ) =0.001 

REFERENCE  SIGNAL  LOOP 
PREDICTION  ERRORS 

ONE  ERROR  ON  BIT  NUMBER  15 
OSCILLATOR 

FREQUENCY  FOR  "0"  BIT  PREDICTION=1010  HERTZ 
FREQUENCY  FOR  "1"  BIT  PREDICTI0N=99O  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCY=2000  HERTZ 
BANDUI DTH= 1 0  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOU  CUTQFF=3  DB 


iyure  V-l.  Simulation  input  data  for  a  case  containing 
a  single  bit  prediction  error. 


1  rom  broadside  and  a  CW  interference  signal  is  incident  at  45°  from 
broadside. 


The  program  output  for  the  case  of  one  bit  prediction  error  is 
shown  in  figure  V-2.  Included  with  the  weight  curves  and  antenna  pat¬ 
tern  is  a  plot  of  the  magnitudes  of  the  desired  and  undesired  parts 
of  the  array  output.  During  the  bit  interval  in  which  the  prediction 
error  occurs,  the  desired  signal  appears  as  interference  to  the  adap¬ 
tive  array.  Therefore,  the  weights  change  quickly  in  an  attempt  to 
null  the  desired  signal.  After  the  bit  interval  containing  the  error 
is  over,  the  weights  begin  to  return  exponentially  back  to  their  steady 
state  values.  We  note  that  tne  array  output  is  amplitude  modulated 
with  an  envelope  similar  to  the  weight  curves  as  a  result  of  the  pre¬ 
diction  error.  The  amount  the  weights  and  the  array  output  are  affect¬ 
ed  by  a  prediction  error  depends  on  the  time  constant  of  the  adaptive 
array.  If  the  array  time  constant  is  several  bit  intervals  long,  the 
array  weights  will  not  change  significantly  during  a  prediction  error. 
Therefore,  with  respect  to  prediction  errors,  the  time  constant  of  the 
adaptive  array  should  be  made  as  long  as  possible.  On  the  other  hand, 
the  array  time  constant  must  be  short  enough  so  that  adaptation  occurs 
in  a  reasonable  period  of  time.  Therefore,  the  choice  of  an  array 
time  constant  is  a  compromise  between  good  tolerance  to  prediction 
errors  and  rapid  array  response. 

In  the  next  set  of  simulation  results,  we  allow  prediction  errors 
to  occur  with  a  probability  of  0.1.  The  CW  interference  is  still 
present  on  the  0  frequency  and  all  other  input  parameters  are  the  same 
as  in  the  previous  case.  A  complete  list  of  the  input  information  is 
given  in  Figure  V-3  and  the  program  output  is  given  in  Figure  V-4. 

We  can  see  in  Figure  V-4  that  bit  prediction  errors  make  the 
weights  less  smooth.  The  slope  of  each  weight  curve  changes  abruptly 
whenever  a  prediction  error  occurs.  From  the  plot  of  the  magnitudes 
of  the  array  output  components,  we  see  that  the  bit  prediction  errors 
tend  to  amplitude  modulate  the  received  signal.  FSK  bit  detection, 
however,  is  relatively  unaffected  by  amplitude  modulation  because  FSK 
bits  are  distinguished  by  frequency.  The  antenna  pattern  taken  at  the 
end  of  the  simulation  period  shows  that  in  spite  of  the  prediction 
errors,  the  array  is  able  to  adapt  properly.  We  see  that  the  adaptive 
array  has  formed  a  null  in  the  direction  of  the  interference. 

The  following  three  sets  of  simulation  results  are  for  prediction 
error  probabilities  of  0.2,  0.3  and  0.4.  Fxcept  for  the  prediction 
error  probabilities,  the  input  parameters  are  unchanged  1 1  o.»  the  pre¬ 
vious  case.  In  Figure  V-5  is  a  list  of  input  data  for  the  case  with 
a  0.2  prediction  error  probability  and  the  corresponding  output  is 
given  in  Figure  V-6.  The  input  data  in  Figure  V-7  and  the  simulation 
results  in  Figure  V- 8  are  for  a  0.3  prediction  error  probability. 

Notice  that  as  the  prediction  error  probability  increases,  the  ampli¬ 
tude  distortion  of  the  array  output  becomes  greater.  Also  we  note 
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adaptive  array  simulation  program  input  data 

INCOMING  SIGNALS] 

DESIRED  F SK  SIGNAL 

"O"  BIT  FREQUENCYS90  HERTZ 
"1*  BIT  FREQUEHCYS010  HERTZ 
P  0  V  £  R  S  D8 

ANGLE-0  DEGREES  FROM  I.  ROADS  IDE 
BIT  R A T E  =  1 0  BITS/SECOND 
MARKOV  TRANSITION  F'ROBAB IL  I  TY=  .  9 

INTERFERENCE  SIGNAL 
FREQUENCY-990  HERTZ 
POU£R= 1 0  DB 

ANGLEM5  DEGREES  FROM  BROADSIDE 
NOISE 

POUER/ELEMENT=-10  DB 

ADAPTIVE  ARRAY 

ELEHENTSS 

SEPARATIONS. 5  DESIRED  SIGNAL  WAVELENGTHS 
CONVERGENCE  FACTORING.  00005 

REFERENCE  SIGNAL  LOOP 
PREDICTION  ERRORS 

PREDICTION  ERROR  PROBABILITY*.! 

OSCILLATOR 

FREQUENCY  FOR  "0"  BIT  PREDICT ION= I  01 0  HERTZ 
FREQUENCY  FOR  "I"  BIT  PREDICTI0N=?90  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCY=2000  HERTZ 
BANDUIDTHMO  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
UNITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOU  CUTOFFS  DB 


Figure  V- 3.  Simulation  input  data  for  a 
prediction  error  probability 
of  0. 1. 
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UDE 


ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 


INCOMING  SIGNALS: 

DESIRED  FSK  SIGNAL 

"0"  BIT  FREQUENCY=9?0  HERTZ 
"1“  BIT  FREQUENCY- 1010  HERTZ 
POUERS  DB 

ANGLE -0  DEGREES  FROM  BROADSIDE 
BIT  RATE  =  1 0  BITS/SECOND 
MARKOV  TRANSITION  PROBAB IL 1 TY- . 8 

INTERFERENCE  SIGNAL 
FREQUENCYS70  HERTZ 
P 0 U E R  = 1 0  DB 

ANGLEM5  DEGREES  FROM  BROADSIDE 
NOISE 

P0UER/ELEMENT=-10  DB 

ADAPTIVE  ARRAY 

ELEMENTS=2 

SEPARATIONS. 5  DESIRED  SIGNAL  WAVELENGTHS 
CONVERGENCE  F ACTOR ( K ) =0 .00005 

REFERENCE  SIGNAL  LOOP 
PREDICTION  ERRORS 

PREDICTION  ERROR  PROBABIL  ITY=. 2 
OSCILLATOR 

FREQUENCY  FOR  "0"  BIT  PREDICTIONSOIO  HERTZ 
FREQUENCY  FOR  "1"  BIT  PREDICT I0N=?90  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCYSOOO  HERTZ 
B ANDU I DTH- 1 0  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOU  CUTOFFS  DB 


riuure  V-5.  Simulation  input  data  for  a  prediction 
error  probability  of  0.2. 


i '  rv'Yv,/’ 

1 1  n  r  |  r  n"r  y~r  rf  t~| 

0  ?0  MO  so 


r  J-  TT  T 

so 


100 


BITS 


-90  -60  -30  0 

DEGREES 


figure  V-6.  Simulation  results  corresponding  to  the 
input  data  of  Figure  V-5. 


ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 


INCOMING  SIGNALS: 

DESIRED  FSK  SIGNAL 

"0“  BIT  FREQUENCY=990  HERTZ 
"1“  BIT  FR£QliENCY=  1010  HERTZ 
P0UER=0  DB 

ANGLE  =  0  DEGREES  FROM  BROADSIDE 
BIT  RATES  0  BITS/SECOND 
MARKOV  TRANSITION  FROBAB  IL I T Y= . 7 

INTERFERENCE  SIGNAL 
FREQUENCY=?90  HERTZ 
POUER= 1 0  DB 

ANGLE-45  DEGREES  FROH  BROADSIDE 
NOISE 

POUER/ELEMENT  =  - 1 0  DB 

ADAPTIVE  ARRAY 

ELEMENTS*? 

SEPARATIONS. 5  DESIRED  SIGNAL  WAVELENGTHS 
CONVERGENCE  FACT  OR ( K ) =0 .00005 

REFERENCE  SIGNAL  LOOP 
PREDICTION  ERRORS 

PREDICTION  ERROR  PROBABILITY*. 3 
OSCILLATOR 

FREQUENCY  FOR  "0"  BIT  PREDICT I0H=1 01 0  HERTZ 
FREQUENCY  FOR  "I"  BIT  PREDICT ION=?90  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCYSOOO  HERTZ 
BANDUI DTH= 1 0  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOLTAGES  DB 
GAIN  BELOW  CUTOFFS  DB 


Fiqure  V-7.  Simulation  input  data  for  the  prediction 
error  pri’chility  uf  0..3. 
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ADAPTIVE  ARRAY  SIMULATION  PROGRAM  INPUT  DATA 

INCOMING  SIGNALS: 

DESIRED  FSK  SIGNAL 

“0“  BIT  FREQUENCYS90  HERTZ 
“1"  BIT  FREQUENCY^ 1010  HERTZ 
POUER-O  DB 

ANGLES  DEGREES  FROM  BROADSIDE 
BIT  RATE-10  BI T S/SECOND 
MARKOV  TRANSITION  PRCBABILI TY= .6 

INTERFERENCE  SIGNAL 
FR'EQUENCY  =  990  HER  IZ 
POUER= 1 0  DB 

ANGLEM5  DEGREES  FROM  BROADSIDE 
NOISE 

P0UER/ELEMENT=-10  DB 

ADAPTIVE  ARRAY 

ELEMENTS=2 

SEPARATIONS. 5  DESIRED  SIGNAL  UAVELENGTHS 
CONVERGENCE  FACT  OR ( K ) =0 . 00005 

REFERENCE  SIGNAL  LOOP 
PREDICTION  ERRORS 

PREDICTION  ERROR  PROBABILITY1 .4 
OSCILLATOR 

FREQUENCY  FOR  "0”  BIT  PREDICT IQN= t 0 1 0  HERTZ 
FREQUENCY  FOR  "I"  BIT  PREDICT I0N=990  HERTZ 
OUTPUT  AMPLITUDES  DB 
LOOP  BANDPASS  FILTER 

CENTER  FREQUENCYSOOO  HERTZ 
BANDUIDTHMO  HERTZ 
CENTER  FREQUENCY  GAINS  DB 
LIMITER 

CUTOFF  VOL T AGES  DB 
GAIN  BELOU  CUTOFFS  DB 


Figure  V-9.  Simulation  input  data  for  a  prediction 
error  probability  of  0.4. 


Simulation  results  corresponding  to  the 
input  data  of  Figure  V-9. 


that  with  the  increase  in  prediction  errors,  the  magnitude  of  the 
interference  at  the  array  output  increases.  In  the  input  data  of 
Figure  V-9,  the  prediction  error  probability  is  set  to  0.4.  The 
corresponding  output,  in  Figure  V-10,  shows  that  the  null  in  the 
antenna  pattern  is  steered  about  halfway  between  the  desired  and 
interference  signals.  With  such  a  high  prediction  error  probability, 
the  array  is  no  longer  able  to  distinguish  the  desired  signal  from 
the  interference.  Therefore,  the  null  in  the  antenna  pattern  is  not 
steered  directly  toward  the  interference. 

Looking  at  the  array  output  curves  in  Figures  V-3  and  Y-5,  we 
see  that  the  average  interference  at  the  array  output  decreases  with 
time  throughout  the  simulation  period.  In  Figures  V-7  and  V-9,  for 
higher  prediction  error  probabi 1 i ties,  the  average  value  of  the  inter¬ 
ference  remains  just  below  the  average  value  of  the  desired  signal. 
Thus  we  conclude  that  for  this  particular  simulation  of  the  system, 
the  prediction  error  should  be  made  0.2  or  less.  This  result,  of 
course,  applies  only  to  the  particular  set  of  system  parameters 
chosen  for  the  simulations. 

The  bit  prediction  error  rate  is  an  important  parameter  in  the 
design  of  the  FSK  system.  As  explained  in  Chapter  III,  the  predic¬ 
tion  error  rate  determines  the  bandwidth  spreading  factor  for  the  FSK 
signal  transmitted  to  the  adaptive  array.  The  prediction  error  pro¬ 
bability  becomes  smaller  as  the  bandwidth  of  the  FSK  signal  increases. 
Therefore,  the  choice  of  the  bit  prediction  error  probability  is  a 
trade-off  between  good  adaptive  array  performance  and  bandwidth  use. 


C.  Reference  Signal  Phase  Error 

Reference  signal  phase  error  can  be  a  problem  in  adaptive  arrays 
that  form  the  reference  signal  from  the  array  output.  By  phase  error, 
we  mean  a  phase  shift  between  the  array  output  and  the  reference  sig¬ 
nal.  The  phase  shift  around  the  reference  signal  loop  is  a  function 
of  frequency.  Therefore,  Doppler  shift  of  the  recei  ed  signal  or 
oscillator  frequency  drift  can  cause  the  reference  signal  to  have  a 
phase  error. 

It  was  shown  by  DiCarlo  and  Compton[13]  that  in  a  one  or  two 
element  adaptive  array,  a  reference  signal  phase  error  causes  the 
weights  to  cycle.  A  subsequent  study  by  D i Ca r 1 o[ 1 4 ]  extended  the 
results  to  include  an  adaptive  array  with  any  number  of  elements  . 

The  study  showed  that  without  interference,  the  cyclic  weights  cai  se 
no  change  in  the  signal  to  noise  ratio  at  the  array  output. 

In  this  report,  we  have  been  studying  the  performance  of  a  two- 
element  FSK  system  in  the  presence  of  CN  interference.  We  are  there¬ 
fore  interested  in  the  effect  of  a  reference  signal  phase  error  on 
the  output  SINR  in  a  two  element  adaptive  array. 
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A  model  of 
f igure  V'll .  We 
In  the  model,  an 
signal  component 


the  reference  signal  phase  shift  problem  is  shown  in 
are  using  the  same  notation  introduced  in  Chapter  II. 
ideal  reference  signal  processor  removes  the  desired 
from  the  array  output  and  produces  a  reference  signal 


r(t)  =  Ar 


sd(t) 

|VdTOT 


(v-i) 


Also  in  the  model  is  a  phase  shifter  which  introduces  a  phase  error, 
ue,  in  the  reference  signal 


r(t) 


=  A 


sd(t) 

r  !sd(t)T 


( V  -  2 ) 


Using  the  equal i ty. 


»  x/u  =  «TXd 


( V-3) 


we  find 

A  X  TW  j  e 
r(t)=  ‘  d  ;;; 

WTXd*)*XdTW 

which  becomes 


r(t)  = 


T  J0e 
Ar  Xd'w  e  e 

y  4>d  w 


(V-4) 


( V  -  5 ) 


The  reference  signal  correlation  vector  is  given  by 

j9e 

A  ♦ ,  W  e  e 

S  =  E{  X*R(t)  )  =  r  r  re:  -  r  ~  .  ( V -6 ) 

VW*T  4>d  W 

The  differential  equation  for  the  LMS  array,  given  by  Equation  (11-45) 
in  Chapter  II,  becomes 
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IDEAL 

REFERENCE  SIGNAL 
PROCESSOR 


ERROR  0e 


Figure  V-ll.  Mathematical  model  used  for  computing 
the  effects  of  reference  signal  loop 
phase  error. 


...  Ar  4- .  W  e 

aw  ,  .  .  a 

+  k<J  w  =  i;  -  -  -  __ 

y  «-!  :du 


(V-7) 


To  get  the  SINR  at  the  array  output,  we  will  find  the  steady  state 
solution  to  the  above  equation  for  the  two  element  case. 

Again  using  notation  introduced  in  Chapter  11 ,  we  calculate 
the  covariance  matrices  for  the  two  element  array. 


.*  =  A  ^ 

'd  Md  U. 


e  u  1 


( V  -  8 ) 


<1i  '  Ai  j«>: 


(V-9) 


y  =  (°?  0 

"  VO  a2 


( V  - 1 0 ) 


The  desired  signal  covariance  matrix,  4^,  has  one  non-zero  eigenvalue 
given  by 


Ad  =  2Ad 


(V-ll) 


The  diagonal  matrix  of  eigenvalues  for  4^  is  defined 


0  2A , 


(V-12) 


We  also  define  a  coordinate  rotation  matrix,  R^,  such  that 
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(V-13) 


R  ,*T  ♦  .  R  .  =  A 
d  d  d  d 


For  the  two  element  adaptive  array,  we  find 


d  ft  1 


( V  - 1 4 ) 


A  new  set  of  rotated  weights, 


(V- i 5) 


are  defined  by  the  equation. 


W  =  R  ,V 
d 


(V- 1 6) 


Substituting  the  above  definition  into  Equation  (V-7)  gives 


...  k  A  <J>  ,R  ,V  e 

R  .  -jy  +k[i>  +<t>  +$  jR  ,V  =  --===■ 

d  dt  d  i  n  d  r. t  +T 

VV*'R  ,  «  .R.V 

a  a  a 


( V  - 1 7 ) 


After  some  manipulation,  we  have 


+T  jop 

dV  kA  R  .  4>,|RjV  ®  *T  *T 

51  =  — L_d_J_d - k  R  *'$  r  v-k  R.  'f.R.V-k*  V  .  (V-18) 

dt  / Tr  Tr — ~ —  d  d  d  did  n  ' 

vV  R  .*  4>  .R  .V 

d  d  d 

Equation  ( V- 1 8 )  can  be  expanded  further  by  using  Equations  ( V- 1 2 ) 
and  (V-13)  into 

I  l*'"*-*^  (I) 


k[o2I+R  *T4-  r  ]  1  !  ( V - 1 9 ) 

0  10  v?| 


Using  Equations  (V-9)  and  ( V - 1 4 ) ,  we  can  expand  the  term 


Rd*T‘iRd 


1  -COS<j> 

js  in<j> 


-jsin<t 

1+COSf 


(V-20) 


where 


t  =  <t 


Vf*i 


(V-21) 


Now,  substituting  Equation  (V-20)  into  Equation  (V- 19),  we  find  a 
pair  of  differential  equations  for  the  rotated  weights 


dv 


Z  2  2 

-  ko  v^-kA.  (l-cos4))v1  +  j kA ^  sin^ 


( V  -  2  2 ) 


dv„  _  v  je  9  9 

dt-  =  /2  “Vr  {vfr  e  2k  Ad  V2  -  k0  v2 


-jkAi?sin4iV1  -  kA^ 2 ( 1  +cos<  )v2  .  (V-23) 


At  this  point,  we  will  assume  a  sinusoidal  steady-state  solu¬ 
tion  for  Equations  (V-2 2)  and  (V-23).  The  assumed  solution  takes 
the  form 


where  V^  and  V2  are  the  magnitudes  of  the  cyclic  weights,  y  is  the 
phase  angle  between  the  two  weights,  and  w  is  the  weight  cycling 
frequency.  c 

Substituting  the  assumed  form  of  the  solution  into  Equations 
(V-22)  and  (V-23),  we  obtain  the  pair  of  equations 


jo.  V1  =  -k(o2+Ai2-Ai2cos^)V1ejY  +  jkA.2sin«.V2  (V-25) 
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jw  V„  --  /2  kA  ,A  e  e  -  k(o2+2A  2+A.2+A.2cos*)V0 
c  .  dr  dii  2 

-j  kAl.2sin<t>VjeJ'Y  .  (V-26) 

Each  equation  above  can  be  divided  into  real  and  imaginary  parts  and 
the  following  four  equations  result: 

i^V^cosy  =  -  k(o2+Ai2-Ai2cosf )V1  sinY  +  kA^sinjV^  (V- 27) 


-<o  V.siny  =  -  k(o2+A^2-A  .2coS(t)V^cosY  (V-28) 

10  V„  =  /2~  k  A  .A  sino„  -  kA  .2sin4>V,cosy  (V-29) 

c  2  a  r  e  i  I 

0  =  /2  kAdArcosoe  -  k(a2+2Ad2+Ai2+Ai2cos<t)V2 

+  kAi2sinfV1 siny  .  (V- 30} 

In  the  following  calculations,  we  will  solve  the  four  equations  above 
by  first  solving  for  the  cycling  frequency,  «c.  Then  we  will  be  able 
to  write  the  other  three  variables,  V],  V?,  and  y,  in  terms  of  uc. 


We  begin  with  Equation  (V-29)  and  obtain  the  relation 
k(o2+A  .2-A  .2cos<>) 

tany  =  -  -1- -  cosy  .  ( V- 31  ) 

(A) 


Now  using  Equation  (V-28)  and  Equation  ( V -31  ) ,  we  find 


koiA^  s i n <}> 

C  OSy  =  ‘  p  p  0  0  "o  ‘ 

W  +k  (o  +Ai  -A.  cos*) 


(V-3 2) 


Substituting  Equation  (V-32)  into  Equation  (V-30),  and  solving  for 
V?  yields 
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( V  -  3  3 ) 


/?  kAdArsinOe[u)c2+k2(o2+Ai2-A.2cos^)2] 

V„  -  '  “  p  p  "p  p  p  '7  7  A  ~  7 

wc[“c  +k  (°  +Ai  -A^  COS*)  +k  a.  sin  *] 

We  substitute  Equations  ( V- 31 )  and  (V- 32)  into  Equation  V-30), 
again  solve  for  V£  obtaining 

/2  kAdArcosee(w2+k2(o2+Ai.2-A.2coS(}>)2) 

V„  =  p  p  p  p  p  p  p  p  p  *  j  p 
[ui  +k  (0  +A-  -A..  cos*)  ] [a  +2Ad  +A^  cos$-kA^  sin  * 

2  2  2 

(a^ai"-A/cos:)]  . 

( V-34 ) 


Equating  Equation  ( V - 33 )  and  Equation  (V-34)  gives  us  an  expression 
for  v 

10  3-io  ktano  (o  +2A  ,  +A .  +A .  cos*) 
c  c  e  dii' 

7  7  7  7  .2  24? 

+  w  (k‘(o  +A.  ^  cos*)  +  k  Ai  sin  *) 

-  k3  tanOe(a2+Ai2-Ai2cos*){(o2+2Ad2+A.2+Ai2cos*) 

(a2+Ai2-A12cos*)  -  A^sin2*}  =  0  .  (V-35) 

Now  by  solving  Equation  (V-35)  above  for  the  cycling  frequency,  we  can 
use  Equation  (V-31 )  to  get  y.  Equation  (V-33)  to  find  V2  and  Equation 
V-32)  to  solve  for  V^ .  The  normalized  weights  then  are  substituted 
into  Equation  (V-16)  to  obtain  the  unrotated  complex  weight  vector  W. 

Fquation  (V-35)  is  a  cubic  with  real  coefficients.  There¬ 
fore  it  will  either  have  one  real  root  and  a  conjugate  pair  of  com¬ 
plex  roots  or  three  real  roots.  In  the  case  of  three  real  roots,  the 
weights  will  oscillate  in  one  of  three  different  frequencies  or  modes 
The  mode  in  which  the  weights  fall  into  depends  upon  the  initial  con¬ 
ditions. 

Equation  (V-35)  was  solved  for  several  values  of  interference 
power  and  the  results  are  shown  in  Figure  V-12.  The  curves  in  the 
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figure  were  drawn  by  computing  the  real  roots  of  the  cubic  equation, 
for  t he  points  where  the  equation  has  three  real  roots,  the  root 
with  the  highest  absolute  value  was  chosen.  At  the  origin,  each  curve 
has  only  one  real  root.  The  case  of  three  real  roots  does  not  occur 
for  any  value  of  phase  error  when  A-j ?  =  5  or  less.  In  the  other  two 
curves  however,  as  the  phase  error  is  increased,  a  point  is  reached 
where  the  cubic  equation  has  three  real  roots.  This  point  appears 
in  the  curves  as  a  slope  discontinuity.  The  discontinuity  appears 
when  three  real  roots  occur  because  one  of  the  three  roots  has  larger 
magnitude  than  the  single  root  at  the  previous  point.  Since  we  are 
plotting  the  largest  root,  a  slope  discontinuity  occurs  in  the  curve. 

We  will  now  study  in  more  detail  a  case  where  the  weights  have 
three  cycling  frequencies.  From  Figure  V - 1 2 ,  we  see  that  for  A,-;=10, 
a  slope  discontinuity  occurs  near  a  phase  error  of  40c .  We  choose 
tiie  point,  f-e=41.65°,  and  calculate  the  three  weight  cycling  fre¬ 
quencies  from  Equation  (V- 35).  For  each  cycling  frequency,  we  can 
find  a  set  of  steady  state  complex  weights  using  Equations  (V-31), 

(V- 3 2)  and  ( V - 33 ) .  A  list  of  the  assumptions,  the  cycling  frequencies, 
and  the  weight  calculations  are  shown  in  Figure  V-13. 

Using  a  fourth-order  Runge-Kutta  a  I gorithm[l 5] ,  we  get  a  numer¬ 
ical  solution  for  the  array  weights  from  Equation  (V-7).  For  the 
case  outlined  in  Figure  V-13,  we  calculate  the  initial  weights,  at 
tO,  to  use  for  starting  points  for  the  Runge-Kutta  algorithm.  Each 
of  the  three  sets  of  initial  conditions  results  in  weights  with  dif¬ 
ferent  cycling  frequencies.  The  weight  curves  for  the  lowest  cycling 
frequency  are  shown  in  Figure  V - 1 4 .  Using  the  initial  conditions 
corresponding  to  the  second  cycling  frequency,  we  obtain  the  results 
of  Figure  V-15.  Finally,  the  solution  for  the  highest  cycling  fre¬ 
quency  is  illustrated  in  Figure  V-16. 

The  relationship  between  the  reference  signal  phase  error  and 
the  array  output  SINR  is  examined  in  Figure  V - 1 7 .  The  curves  were 
drawn  by  first  computing  the  cycl ing  frequencies  and  weights  for  each 
value  of  phase  error.  Then  the  components  of  the  array  output  due 
to  desired  signal,  interference,  and  noise  were  calculated.  For  the 
cases  where  three  cyclinq  frequencies  occurred,  the  freauenev  resultina 
in  min.imunj  SINR  was  chosen.  In  Figure  V-17,  we  see  the  array  output 
SINR  in  decibels  for  several  values  of  interference  power.  For  the 
case  of  no  interference,  we  see  that  the  SINR  remains  constant  with 
increasing  phase  error.  When  interference  is  added,  however,  an  in- 
<  rease  in  phase  error  causes  a  reduction  in  SINR.  As  more  interference 
power  is  added,  the  reduction  in  SINR  becomes  more  pronounced.  In 
the  last  two  curves  of  Figure  V-17,  there  is  an  almost  vertical  drop 
in  SINR  as  the  phase  error  is  increased  past  a  certain  point.  This 
SINR  drop  occurs  at  the  point  where  the  solution  to  the  equation  for 
the  weight  cycling  frequency  goes  from  one  real  root  to  three  real 
roots.  These  curves  show  that  for  high  interference  levels,  a  large 
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Mode  1 . 


=  2.865  rad/sec 


Figure 


^C] 

V1  =  0.1554 
V2  =  0.1116 
y  =  59.903° 


J  -Cl  ^ 

W,  --  (0.09301-jO. 09508)  e 

j:.dt 

W2  =  (0.0238-jO. 09508)  e 


Mode  2. 


C2  =  5.560  rad/sec 


V]  =  0.08433 
V2  =  0.07885 
Y  =  41.460° 


W]  -  (0.01 119-jO. 03962)  e 
W 2  =  (-0.10032-jO. 03962)  e 


Mode  3.  =  8.527  rad/sec 

V]  =  0.05388 
V2  =  0.06674 
Y  =  30.107° 

j 

W1  =  (0.01423-jO. 01911)  e 
W ?  -  ( -0.08015-jO. 01911)  o 


V- 1 3.  Data  and  results  for  a  three-mode  example. 
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reference  signal  phase  error  can  cause  a  serious  reduction  in  the  per¬ 
formance  of  the  adaptive  array. 

We  have  seen,  in  this  section  the  effects  of  a  reference  signal 
phase  error  in  the  adaptive  array.  First,  we  showed  that  a  phase 
error  causes  the  steady  state  array  weights  to  cycle.  Under  some 
conditions  there  is  only  one  cycling  frequency  but  sometimes  Equation 
(V-35)  gives  three  real  cycling  frequencies.  When  no  interference  is 
present,  we  showed  that  the  phase  error  has  no  effect  on  SINR.  How¬ 
ever,  we  showed  that  in  the  presence  of  interference,  the  phase  error 
a  significant  reduction  in  SINR.  Therefore  we  conclude  that  in  de¬ 
signing  the  reference  signal  loop,  it  is  important  to  keep  the  phase 
error  as  small  as  possible. 


D.  Signal  Acquisition  and  Synchronization 

We  assumed  in  our  discussion  of  the  reference  signal  loop  in 
Chapter  IV  that  bit  predictions  were  available  during  the  array  lock¬ 
up  period.  We  also  assumed  that  the  bit  predictions  were  synchron¬ 
ized  with  the  incoming  desired  signal  bits.  These  two  assumptions 
imply  that  we  can  detect  the  desired  signal  and  extract  timing  in¬ 
formation  during  the  array  lock-up  period.  In  this  section  we  will 
show  how  the  array  responds  in  the  early  part  of  the  acquisition 
period,  before  the  bit  predictions  and  the  reference  signal  are  avail¬ 
able.  It  will  be  shown  that  even  without  the  reference  signal  present, 
the  array  maintains  a  good  SINR  during  the  initial  lock-up  period. 

The  reference  signal  loop  requires  properly  timed  bit  predic¬ 
tions  to  produce  a  reference  signal  that  correlates  with  the  desired 
signal.  During  the  early  part  of  lock-up,  the  interference  power  at 
the  array  output  may  exceed  the  desired  signal  power.  Therefore,  it 
may  not  be  possible  to  detect  the  desired  signal  bits  and  make  bit 
predictions.  When  the  bit  predictions  are  not  available,  we  will 
assume  that  the  output  of  the  reference  signal  loop  does  not  corre¬ 
late  with  the  desired  signal.  Therefore  the  reference  signal  corre¬ 
lation  vector  becomes 


S  =  EIX(t)r(t)l  =  0 


( V  -  36 ) 


When  the  reference  signal  does  not  correlate  with  the  desired  signal , 
the  differential  equation  for  the  LMS  array  is 


dW 

dt 


k  iW 


( V  -  37  ) 
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Figure  V - 1 8  shows  curves  of  SINR,  output  desired  signal  power 
(Ptf),  and  output  interference  power  (Pj)  for  the  period  where  no  bit 
predictions  are  available.  The  curves  were  obtained  by  using  a  Runge- 
Kutta  algorithm  to  solve  Equation  ( V - 37 )  for  the  array  weights.  The 
weights  were  then  used  to  compute  the  signal  powers  at  the  array  out¬ 
put. 

We  see  in  Figure  V- 1 8  that  the  adaptive  array  does  power  inver¬ 
sion  when  no  reference  signal  is  present.  By  power  inversion,  we 
mean  that  after  the  initial  transient, 


Therefore,  when  the  incident  interference  power  is  much  greater  than 
the  incident  desired  signal  power,  the  array  will  quickly  give  a 
good  S 1  NR  at  the  array  output.  The  absolute  signal  levels  however 
decrease  with  time.  Therefore  acquisition  of  the  signal  and  timing 
must  take  place  after  a  good  SINR  is  established  and  before  the 
signals  become  too  small  to  be  detected. 

The  FSK  detector  used  in  the  system  described  here  is  asso-ed  to 
be  noncoherent.  The  problem  of  bit  synchronization  for  noncoherent 
has  been  studied  by  Resso  [16].  We  will  not  go  into  the  problem  of 
bit  timing  in  this  report.  The  reader  is  referred  to  Resso  for  details. 

We  have  seen  that  in  the  presence  of  large  interference  levels, 
the  power  inversion  property  of  the  LMS  array  gives  us  a  time  period 
in  which  to  acquire  the  desired  signal  and  timing.  We  will  assume 
that  this  time  pe'iod  is  long  enough  that  bit  detection  and  synchroni¬ 
zation  will  not  be  a  serious  problem. 
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CHAPTER  IV 
CONCLUSION 


We  have  examined  the  use  of  an  adaptive  array  in  an  FSK  commu¬ 
nications  system.  The  adaptive  array  is  used  to  protect  the  commu¬ 
nications  system  from  interference.  Interference  protection  is 
especially  valuable  in  military  communications  for  example,  where 
deliberate  jamming  may  be  a  factor. 

The  FSK  system  presented  in  this  report  uses  an  LMS  adaptive 
array.  As  explained  in  Chapter  II,  the  LMS  array  minimizes  the  mean 
square  error  between  the  array  output  and  a  reference  signal .  There¬ 
fore,  to  use  the  LMS  array  in  an  FSK  system,  we  must  generate  an  FSK 
reference  signal . 

We  described  a  method  for  forming  an  FSK  reference  signal  based 
on  the  use  of  bit  predictions.  In  Chapter  III  we  showed  that  by 
making  the  desired  bit  stream  first  order  Markov,  bit  predictions 
can  be  made  at  the  FSK  receiver. 

In  Chapter  IV,  we  described  a  reference  signal  loop  which  uses 
the  array  output  and  the  bit  predictions  to  form  a  reference  signal. 
When  the  reference  signal  loop  operates  properly,  the  desired  signal 
component  of  the  array  output  is  fed  back  to  the  adaptive  array  as 
the  reference  signal . 

Several  possible  problems  associated  with  the  FSK  system  were 
discussed  in  Chapter  IV.  We  first  looked  at  the  effects  of  bit  pre¬ 
diction  errors.  Then  we  considered  the  problems  of  reference  signal 
phase  error  and  the  resulting  cyclic  weights.  Finally  we  looked  at 
signal  acquisition  and  bit  synchronization  during  lock-up. 

In  conclusion,  we  have  presented  an  adaptive  array  system  for 
use  in  FSK  communications.  The  system  operation  was  demonstrated 
with  the  use  of  a  computer  simulation  program.  We  found  that  within 
the  limits  discussed  in  this  report,  the  FSK  system  developed  here 
is  a  viable  system  for  use  in  FSK  communications. 
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APPENDIX 


0010  c 

0020  C 
0030  C 
0040  C 
0050  C 
0060  C 
0070  C 
0080  C 
0090  C 
0100  C 
0110  C 
0120  C 
0130 
0140 
0150 
0160 
0170 

0180  C 
0190  C 
0200  C 
0210  C 
0220 
0230 
0240 
0250 

0260  C 

0270 

0280 

0290 

0300  C 
0310 

0320  C 
0330 

0340  C 
0350 
0360 
0370 

0380  C 

0390 

0400 


SIMULATION  PROGRAM  FOR  AN  ADAPTIVE  ARRAY  SYSTEM  USING 
FSK  MODULATION. 

THIS  PROGRAM  SIMULATES  A  TWO  ELEMENT  LMS  ADAPTIVE  ARRAY, 

A  MARKOV  SOURCE  FOR  THE  DESIRED  SIGNAL  BITS,  A  BIT  PRED¬ 
ICTOR  CIRCUIT,  AND  A  REFERENCE  SIGNAL  LOOP. 

THE  SIMULATION  IS  DONE  IN  DISCRETE  TIME  BY  SAMPLING  EACH 
SIGNAL  AT  A  RATE  EQUAL  TO  FIVE  TIMES  THE  CENTER  FREQUENCY 
OF  THE  DESIRED  FSK  SIGNAL. 

THE  PROGRAM  CALCULATES  AND  STORES  THE  WEIGHT  SAMPLES  AND 
THE  OUTPUT  CONSISTS  OF  PLOTS  OF  THE  ARRRY  WEIGHTS  AND  AN 
ANTENNA  PATTERN  CALCULATED  AT  THE  FINAL  TIME  SAMPLE. 

DIMENSION  T A < 500 ) , U1  ( 500) ,U2 < 500 ) , U3( 500) , U4 ( 500 ) , T API 500 ) 
DIMENSION  AS0UT(500),AI0UT (500) 

DIMENSION  AO T < 5 )  ,SN(5)  ,S0(5) 

DIMENSION  ANG(361),DBS(36U 
PI =3 . 14159265359 


INPUT  DATA 

DESIRED  SIGNAL  PARAMETERS 
FS0=990. 

FS1=1010. 

AS=1 .0 
THS=0. 

CU  INTERFERENCE  PARAMETERS 
AI  =  1  . 

FI=990. 

THI=45. 

MARKOV  BIT  CHANGE  PROBABILITY 
PA=1 . 

SAMPLING  RATE 
FSAM=5000. 

OSCILLATOR  PARAMETERS 
F0S0=1 010. 

F0S1 =990. 

DOSsO. 

LIMITER  PARAMETERS 
GL*2. 

C0FF=1 . 
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0410 

C 

REFERENCE  LOOP  BANDPASS  FILTER  PARAMETERS 

0420 

FC=2000. 

0430 

8U  =  5. 

0440 

C 

LMS  ARRAY  FEEDBACK  LOOPS  GAIN  CONSTANT 

0450 

C0NV=0. 00005 

0460 

C 

ANTENNA  SEPARATION  IN  DESIRED  SIGNAL  UAVELENGTHS 

0470 

AL=0.5 

0480 

c 

THERMAL  NOISE  POUER 

0490 

S2= .  1 

0500 

c 

0510 

c 

INITIAL  VALUES 

0520 

c 

0530 

c 

SEED  FOR  THE  RANDOM  NUMBER  ROUTINE 

0540 

1X-4021 947 

0550 

c 

MEAN  VALUE  OF  THE  THERMAL  NOISE 

0560 

AM=0.0 

0570 

c 

INITIAL  TRANSMITTED  BIT 

0580 

rt  S I G  =  1 

0590 

c 

INITIAL  BIT  PREDICTION 

0600 

MOSC= 1 

0610 

c 

INITIAL  i-IEIGNT  VALUES 

0620 

Ul<1)=1.0 

0630 

U2( 1  )* 1 .0 

0640 

U3( 1 )-1 .0 

0650 

U4(1 )  =  1  .0 

0660 

c 

INITIAL  TIME  VALUE 

0670 

TA(1 1=0. 

0680 

c 

INITIAL  VALUES  FOR  TRIG  FUNCTIONS.  TRIG  VALUES  ARE 

0690 

c 

TO  PREVENT  ERROR  CAUSED  BY  COMPUTING  TRIG  FUNCTIONS 

0700 

c 

LARGE  VALUES. 

0710 

0SIN0=0. 

0720 

OCOSO  =  1 . 

0730 

0SIN1 =0. 

0740 

OCOS1=1. 

0750 

DSIN0=0. 

0760 

DC0S0=1. 

0770 

DSIN1=0. 

0780 

DC0S1 =1 . 

0790 

SINUI=0. 

0800 

C0SUI=1. 

0810 

c 

UNIT  CONVERSIONS 

0820 

USO=2.*PI*FSO 

0830 

US1 =2.*PI*FS1 

0840 

UI=2.*PI*FI 

0850 

DT=1 ./FSAM 

0860 

UOSO=2.*PI*FOSO 

0870 

,  U0S1  =»2.*PI*F0S1 

0880 

AOS=PI*DOS71 80. 

0890 

MC*2.*PI+FC 

0900 

UBU=2.*PI*BW 
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0910 

0  v  7  o 

0930 
0940 
0950 
0960 
0970 
0980 
O990 
1000 
10 '0 
; A 
1  030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 

I  t  30 

I I  40 
50 

1160 
tl  70 
1180 
1190 
1200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 


C  INITIALIZE  STACK  VECTORS 
DO  10  1=1,5 
AOT  < I ) =0. 

SN( I ) =0. 

10  S0(I)=0. 

T=0 .0 
IUT  =  1 

C  INITIALIZE  SCa  vEIGHTS 
UP1=U1<1) 

UP2=U2  < 1 ) 

UP  3  =  U  3  < 1 ) 

UP  4  - U 4  <  1  ) 

C  SIGNAL  SPACE  ANGLES 
U5C= ( US0+US1  )/2. 

BTSO=2.*PI+AL*$IN<THS*PI/180.)*USO/USC 
BTS1=2.*PI*AL*SIN(THS+P 1/180.) *US1 /USC 
BTSA= ( BTS0+BTS1 )/2. 
BTI=2.*PI*AL*SIN<THI*PI/180.)*WI/USC 
C  BIT  COUNTER 

DO  40  «=1 , 1 00 

C  PREDICT  OFPOSITE  OF  LAST  BIT 
IF < NSIG. EQ. 1 )  MQSC=0 
IF(HSIG.EQ.O)  N0SC=1 

C  MARKOV  SOURCE;  FIRST  ORDER  SYMMETRIC 
CALL  RANDU<IX,IY,RNUM) 

IX=IY 

IF  (RNUM.LT  .F'A.AND.MSIG.EQ.O)  MSIGP=1 
IFIRNUN.LT.PA.AND.HSIG.EQ.I  )  MSIGP=0 
MSIG=MSIGP 

C  RANDOM  PHASE  FOR  INTERFERENCE  BIT 
CALL  RANDU< IX, I Y,RNP) 

I X  =  I Y 

RNP=RNP*2.*PI 

C  UEIGHTS  SAVED  PER  BIT  COUNTER 
DO  35  1=1,5 

C  1TTERATI ONS  PER  WEIGHT  SAVED  COUNTER 
DO  30  J=1,100 

C  DESIRED  SIGNAL  CALCULATIONS 

DSIN0P  =  DS1N0+C0S(WSQ  +  DT )  +  DCGSO*SIN (USO +DT  > 
BCOSOP=DCOSO*COS(USO«DT)-DSINO«SIN(USO+DT) 
DSIN1P=DSIN1*C0S<US!*DT)tDC0S1*SIN(US1*DT) 
DC0S1P=DC0S1*C0S(US1*DT)-DSIN1*SIN(US1+DT) 

IFIMSIG.EQ. 1 )  GO  TO  20 
XS1 = AS»DCOSOP 
XS2=AS*DSIN0P 

XS3  =  AS,M  DCOSO+COSI WSO*DT-BTSO) -DSINO+SIN ( USO • DT  BT  SO > ) 
XS4=AS* ( BS1NO*COS(USO*DT-BTSO)+DCOSO+SIN(USO*DT-BTSO) > 
GO  TO  13 

20  XS1=AS*BC0S1P 
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1 4  00 
MIC 
14  20 

14  30 
t  -MO 
I  4  50 
1460 
1470 
1480 
1  4  VO 
1500 
1510 
1520 

15  30 
1540 
15  50 
1540 
1570 
1580 
1590 
1600 
1610 
1420 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
17  90 
1800 
1810 
1820 
1830 
1840 
1850 
1860 
1870 
1880 
1890 


XS2= AS*DSIH1 P 

XS3=AS*(I>C0SHC0S(US1*IlT-BTSn-DSINt+SIN(US1*IiT-»TS1 ) ) 
XS4-AS*(DSIN1*C0S(US1*DT-BTS1 ) +DC0S1 *SIN( WS1 *IlT-BTSl ) ) 

13  DSINO=DSINOP 
DCOSOrDCOSOP 
DSI«1  =  DSINtP 
DC  OS  1 =DC0S1P 

C  INTERfER£NCE  SIGNAL  CALCULATIONS 

XI1=AI*(C0SUI*C0S<UI *DT+RNP)-S1NUI*SIN(UI*DT+RNP> ) 
XI2^AI*(SINUI*COS(UI*DT+RNP)+COSUI*SIN(UIrDT+RNP)) 
XI3=Al4(C0SUI*C0S(UI*DT-BTI+RNP)-SINUI*SIN(UI+DT'-bTI+RNP) ) 
XI4=AI*<SINUI*C0S(UI*DT-BT1+KNP)+C0SUI*SIN(UI *DT -BT I ♦RNP ) ) 


SINUIP-5INU1*C0S(UI*BT)+CCSL'I*SIN(UI*DT) 

COSUIP-COSUI*COS(U2»DT)-S'INOI»SIN(UI*IiT) 

SINUI=S1NUIP 

C0SUI=C0SU2P  : 

C  RANDOM  NOISE  jj 

CALL  GAUSS1  IX,S2,AM,PN0S1  )  \\ 

CALL  GAUSS  ( I X  , S2 , AM ,RN0S2  )  |! 

CALL  GAUSS<IX,S2,AM,RN0S3)  i! 

CALL  GAUSS(IX,S2,AM,RN0S4)  jj 

C  COMPOSITE  SIGNALS  I] 

XI  =XS  1  +X 1 1  +RN0S1  !; 

X2-XS2+XI2+RN0S2  j 

X3=XS3+XI3+RN0S3 

X4=XS4+XI4+RN0S4  j 

C  ARRAY  OUTPUT 

ASIG-UP 1  ■*xi+UP2*X2+UP3-fX3  +  UP4-*X4  !' 

UFACS=SQRT(UPI»42+UP2**2+UP3*»2+UP4i'+2+2.*(UP1*UP3+UP2tUP4) 

&BTSA)+2.*(UP2*UP3-UP1*UF'4)*SIN(BTSA)) 

UFACI=SQRT(UP1**2+UP2**2+UP3«+2+UP4*»2+2.*<UP1*UP3+UF'2*UP4) 

XBTI)»2.*(UP2*UP3-UP1*UP4)*SIN<BTI)) 

THUI=ATAN2(UP2+UP3*SIN(BTI)+UP4+C0S(BTI>,  ! 

iUP1+UP3*COS(BTI)-UP4*SlN(BTI>)  j 

THUO=ATAN2(UP2+UP3*SIN(BTSO)+UP4*COS(BTSO), 
SUF1+UP3*COS<BTSO)-UP4*S1N(BTSO>> 

THU1=ATAN2(UP2+UP3*SIN(BTS1 ) +UP4»C0S ( BTSU , 

SUF'1+UP3*C0S(BTS1  )-UP4*SIN(BTS1 )) 

C  BEGIN  CALCULATIONS  FOR  LIMITER  OUTPUT.  THE  FOLLOUING  CODE 

C  ALLOUS  FOR  WEAKER  SIGNAL  ATTENUATION  IN  THE  LIMITER. 

C  REFERENCE  COMPONENT  DUE  TO  INTERFERENCE 

AIP=AI«UFACI 
ULI=UI+UOSO 

IF (MOSC . EQ. 1 )  ULI-UI*U0S1 
CALL  LFILT(AIFAC,ULI ,UC,UBU) 

AIP=AIP*AIFAC+GL 

C  REFERENCE  SIGNAL  COMPONENT  DUE  TO  SIGNAL 

ASP=AS*UFACS 
IF<MSIG.E0. 1 )  GO  TO  61 
ULSrUSO+UOSO 
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1<?00 

IF (N0SC.E9. 1 )  ULS=gS0*U0S1 

1910 

CALL  LFILT (ASFAC,ULSfUC,U8W> 

1920 

ASP=ASP*ASFAC*GL 

1930 

GO  TO  62 

1940 

61 

ULSrUS1 +U0S0 

1950 

IF (HOSC.EG. 1 )  ULS=US1 +UQS1 

I960 

CALL  LFILT1 ASFAC,ULS,UC,UBU) 

1970 

ASP=ASP4ASFAC*GL 

1980 

42 

CONTINUE 

19<>0 

CQFF2=2.+CGFF 

2000 

IFIASP.LE .C0FF.Ai9D.AIP.LE.C0FF)  GO  TO  78 

2010 

IF i ASP.LE.CQFF . AND. AIP.GE.C0FF2)  GO  TO  63 

20  20 

GO  TO  64 

20?" 

63 

ASF=ASP/2. 

20  4  C 

A1P=C0FF 

2050 

GO  TO  78 

2060 

64 

CONTINUE 

2070 

IF (AIP.LE.COFF .ANP. ASP.6E.C0FF2)  GO  TO  65 

O 

no 

o 

GO  TO  66 

2090 

65 

ASP=C0FF 

2100 

AIP-AIP/2. 

2110 

GO  TO  78 

2120 

66 

CONTINUE 

2130 

IF(AIP.LE.C0FF.ANP.ASP.GT.CQFF.AND.ASP.LT.C0FF2> 

2140 

GO  TO  68 

2150 

67 

ASP^COFF 

2160 

AIP=COFF*AIP/ASP 

21  70 

GO  TO  78 

2 1  BO 

68 

CONTINUE 

ziyo 

IF ( ASP .LE.COFF .AND. A IP .GT .COFF .AND. A IP. IT .C0FF2) 

2200 

GO  TO  70 

2210 

69 

ASP'COh  F  *ASP/AIP 

2220 

AIP=C0FF 

2230 

GO  TO  78 

2240 

70 

CONTINUE 

2250 

RAT I0=AIP/ASP 

2260 

IF(RATIQ.LT.0.5)  GO  TO  71 

2270 

GO  TO  72 

2280 

71 

ASP=CQFF 

2290 

AIP=C0FF/2. 

2300 

GO  TO  78 

2310 

72 

CONTINUE 

’320 

IFIRATIO.GE. 0.5. AND. RATIO. LT.I.)  GO  TO  73 

* » o 

GO  TO  74 

J  40 

73 

ASP-COFF 

.'350 

AIP=COFF+RATIO 

<60 

GO  TO  78 

i  0 

74 

CONTINUE 

■  “O 

IF (RATIO. GE.1. ANO.RATIO.LT. 2.)  GO  TO  75 

GO  TO  74 


78 


2400 

75 

AS  P= C  0  F F /RAT  10 

.'410 

AIP=CQFF 

2420 

GO  TO  78 

24  30 

76 

CONTINUE 

2440 

IF (RATIO. GE .2. )  GO  TO  77 

2450 

GO  TO  78 

2460 

77 

ASP=CQFF/2. 

2470 

A 1 P= COF  F 

2480 

78 

CONTINUE 

2490 

REFI=AIP+(CQSUI*COS (RNP-THUI )-SINUI*SIN (RNP-THUI > ! 

2500 

IF (M51G.EQ. 1 )  GO  TO  79 

2510 

REFS=ASP*<DCOSO*COS(THUO)+USINO*SIN(THUO) ) 

252  0 

GO  TO  80 

:530 

79 

REFS=ASP*(PCQS1*C0Si THU  1) +HS INI *S INC  THU  1 ) > 

5  4  rj 

80 

CONTINUE 

.550 

C 

TOTAL  REFERENCE  SIGNAL 

2560 

REF=REFS+REFI 

2570 

C  ERROR  SIGNAL  VALUE 

2580 

ER=ASIG-REF 

2590 

C  CALCULATE  NEU  WEIGHTS 

2600 

UP1=UP1-2.*C0NV*ER*Xt 

2610 

UP2=UP2-2.*C0NV*ER*X2 

2620 

UP3=UP3-2. +CQNV*ER*X3 

2630 

UF4:UP4-2.*C0NV*ER*X4 

2640 

T=T+DT 

2650 

30 

CONTINUE 

2660 

C  SAVE  EVERT  20TH  WEIGHT  VALUE  FOR  PLOT 

2670 

IUT=IUT+1 

2680 

IF ( IUT -GT .500)  GO  TO  35 

2690 

U1 ( IUT ) =  UP  1 

2700 

U2 ( IUT ) ~UP2 

2710 

U3( I U  T )  =  UP  3 

2770 

U4 ( IUT )=UP4 

27j0 

ASOUT ( IUT)=AS*UFACS 

2740 

AI0WT(IUT)=A1*UFACI 

2750 

TA ( IUT  )  =  T-DT 

2760 

35 

CONTINUE 

2770 

40 

CONTINUE 

2780 

DO  99  1=1,500 

2790 

99 

TAP(I)=TA<I)*10. 

2800 

C 

PLOT  RESULTS 

2810 

C 

GENERAL  PLOTTING  ROUTINE  USED 

2820 

CALL  PLTPKM < TAP, U1, 500, 500, -1,0,1 ,.25 > 

2830 

CALL  PLTPKMI TAP, U2, 500, 500,-1 ,0,1 ,.25) 

2840 

CALL  PLTPKMI TAP,U3, 500, 500,-1 ,(,1,.25> 

2850 

CALL  PLTPKMI  TAP,  U4, 500, 500,-1  ,-.’,1  ,.25) 

2860 

CALL  PLTPKMI TAP, ASOUT, 500, 500, -1 ,0,2, .25) 

2870 

CALL  PLTPKMI TAP, A I OUT ,500 ,500, -2,0, 2 , .25) 

2860 

CALL  PATPLT(AL,U1(500),U2(500) ,U3 (500) ,U4 ( 500) , ANG 

2890 

CALL  EXIT 

2900 

END 

79 


i 


2910 

C 

2920 

c 

SUBROUTINE  FOR  ANTENNA  PATTERN  PLOTS 

2930 

c 

2940 

SUBROUTINE  PATPLT(AL,U1 ,U2,U3,U4,ANG,BBS,N> 

2950 

DIMENSION  ANG( N) , DBS(N) 

2960 

PI -3. 1 4159265359 

2970 

NP=N-1 

2980 

DO  5  1  =  1,  N 

2990 

AN6< I)=PI* (FLOAT (1-1 >/NP~Q.5> 

3000 

RETA=2.*PI*AL+SIN(AN6(I>) 

3010 

5 

DBS ( I ) =U1  +U1 +U2*U2+U3*U3+U4*U4+2.* (W1 +W3+U2+U4 1* 

3020 

SCGS<&ETA>+2.*(U2*U3-Ul*U4>*SIN<BETA> 

3030 

AM  A  X  =  0 . 

3040 

DO  10  1  =  1, N 

3050 

IFUiBS(I).LT.AHAX)  60  TO  10 

3060 

ArtAX--DBS(I) 

3070 

10 

CONTINUE 

3080 

DO  15  1=1, N 

3090 

DBS  1 1 )  =  1 0 . + ALQG1 01 ABS IBBST 1 1  / ANAX+ 1  -OE-30 >  > 

3100 

15 

ANG(I)= AN 6(11*180. /PI 

3110 

(1N=-N 

3120 

CALL  PLTPKM(ANG,DBS,N,KN,-1 ,0,1 ,.25) 

3130 

RETURN 

3140 

END 

3150 

C 

3160 

C 

SUBROUTINE  FOR  GAUSSIAN  RANDOM  NOISE 

3170 

C 

3180 

SUBROUTINE  GAUSS(IX,S2,AH,V2) 

3190 

A=0.0 

3200 

DO  50  1=1,12 

3210 

CALL  RANDU ( IX, I Y, V ) 

3220 

IX=IY 

3230 

50 

A=A+Y 

3240 

V2=<A-6.0)*S2+AN 

3250 

RETURN 

3260 

END 

3270 

C 

3260 

C 

FUNCTION  FOR  RANDOH  NUMBER  GENERATION 

3290 

C 

3300 

SUBROUTINE  RANDUUX) 

3310 

I Y=IX*6553? 

3320 

1F(1Y)5,6,6 

3330 

5 

IY=IY+21 47483647+1 

3340 

6 

YFL=IY 

3350 

YFL=YFL* .465661 3E'9 

3360 

RETURN 

33  ’0 

END 
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33NO  ( 

3390  ! 

2400  I 

34)0 

1420 

3430 

3440 

34‘jO 

3460 

34  '0 

3480 

3490 

3S00 

3  j  1 0 

3‘ 20 


LOOP  BANDPASS  FILTER  SUBROUTINE 

SUBROUTINE  LFILT(FAC,UL,UOfUBU) 

COMPLEX  G,H,S 

AL=0.5t>UBU 

S=CMPLX(0.,U0) 

GrS/  ( S*S+2. *AL*S+ AL*AL+UO*UO) 

G=t . /G 
AK-CABS(C) 

S=CMPLX(0. ,UL) 

H=S*AK/(S4S+2.*AL*S+AL*AL+U0«U0) 
F  AC  =  CABS ( H ) 

RETURN 

END 
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